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ABSTRACT 


A procedural approach to the design and implementation 
of a computer simulation model is presented, with a simula- 
tion model description and computer code. The discussion 
and example are intended as reference material for the 
computer science and computer war gaming courses offered 
at the United States Naval Postgraduate School. The 
Sample computer simulation was designed for the statisti- 
cal analysis of the comparative effectiveness of different 
ASW helicopter search tactics in a variety of tactical and 
physical environments. This simulation has available a 
wide range of input parameters and is applicable to all 
ASW] helicopters and any search plan employing ten helicop-~ 
ters or less. The accompanying FORTRAN computer code was 
written for the CDC 1604 computer, and is adaptable to the 
IBM 7090-94 by the inclusion of the appropriate control 


cards and random number generator. 
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CHAPTER I 


INTRODUCTION 


m1 Objective. 
The objective of this thesis is to offer the beginning 


student of computer war game simulation an easily understood 
example of what a computer simulation is and how one might 
go about constructing one. The first year Operations 
Analysis students of the U. S. Naval Postgraduate School 
constitute the particular audience for which the thesis is 
intended. Accordingly, no more than a basic understanding 
of probability and the Monte Carlo Method is presupposed. 
Fundamentals of War Gaming by F. J» McHugh [1] and A Survey 
of Historical Developments in War Gaming by John Young [2] 
are recommended to the reader who has not studied manual 
war gaming methods or is unfamiliar with the historical 
development of war games. 

The particular method chosen by the author to reach 
the stated objective consists of a sample computer simula= 
tion and a proposed plan for designing and implementing 2 
computer simulation model. The proposed simulation design 
methodology is essentially one solution to the major problems 
confronting a novice computer simulation designer. 

The documentation and computer code for a helicopter 
anti-submarine warfare search simulation are included in Chapter 
VI and Appendix A respectively. This simulation, subsequently 
referred to as AHS-l1, is cited as an example throughout this 


thesis. AHS-1 was written as an example because 
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it treats an operation of general interest in the U. S. Navy, 
ives; @ Seurch Tor an evading sudmarine; and is of limited 
enough scope to be easily understood. This particular subd- 
ject was selected as a result of the author's past experience 
as an Anti-Submarine Warfare Officer in helicopter squadrons. 
AHS-1 is intended as an analytic tool for the evaluation of 
the relative effectiveness of different helicopter search 
plans in similar environmental and tactical circumstances, 

It is hoped that units within the U. S. Navy, whose responsie- 
bility it is to conduct such analysis, will consider using 
the simulation; either as written or with any necessary modi- 
fications. With this application in mind it should be noted 
that a hypothetical distribution of helicopter figure of merit 
was incorporated in the simulation. This allowed the thesis 
to be unclassified. Rewriting this particular portion of the 
coding is a relatively minor modification. 

1.2 Definitions. 

Some confusion exists throughout the literature in the 
definition of terms used in War gaming and simulation in gene 
eral. This is perhaps partly attributable to the diverse 
backgrounds of the individuals engaged in the field. War 
gamers include people from nearly every branch of science and 
engineering as well as professional military men. This con- 
fusion is possibly also due in part to a tendency to attach 
special meaning to certain terms according to the particular 
activity in which the individual is engazed. The following 
definitions will be adhered to in the remainder of this thesis 


and were chosen for simplicity as well 1s renerality. 
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War Game has been applied to activities ranging in com- 
plexity from two opponents at a sand table to strike exer- 
cises involvins thousands of men and millions of dollars 
worth of equipment. A. We Pennington chooses to leave the 
term undefined because of the almost universal disasreement 
on a definition. 3 | F. do McHugh defines a War game as 
"a simulation, in accordance with predetermined rules, data, 
and procedures, of selected aspects of a conflict situation." 
[a] This definition is favored by the author because it is 
clear and concise and does not tend to imply any special type 
of war game, size or mode of play. The words "conflict sit- 
uation" should be particularly noted. This is the primary 
distinction between a@ game and other types of simulation. In 
accordance with current usage, war game will not include ex- 
ercises in which real forces or equipment are employed. 

Because of its size and complexity, the real world, or 
even selected aspects of the real world, is usually impossible 
to study in detail in the laboratory. Consequently, scientists 


must usually resort to the study of a model of a real life 





Situation. This model may be iconic, such as a model airfoil 
constructed by the aeronautical engineer for study in a wind 
tunnel, or it may be symbolic. A stick and eumdrop model of 

an atom is a symbolic model, as is the familiar equation F = ma. 
The latter is an example of a mathematical model and together 
with several other equations, or models, is the model for 
Newtonian mechanics. In essence, a model is a representation 
of some aspect of real life on which experiments can be pere- 


formed; and from the behavior of the model certain laws or 
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general rules governing the real world situation may be ine 
ferred. It is characteristic of every model thst the result= 
ing conclusions adout the real world will never be any more 
valid than the model itself. This thesis will be concerned 
with mathematical models. 

Mathematical models can be further classified as deter~ 
ministic or stochastic. The equation F = ma is a determin- 


istic model in that if mass and acceleration are ‘cnowm tne 





force is uniquely determined. If, however, the measurement 
of acceleration is subject to certain random errors the above 
equation could be written as 
F = m(ate) 

Where a is the measured value of acceleration and e re- 
presents a random measurement error. In this example, F is 
@ random variable, and the equation is a stochastic model 
representing an abstract concept which is called force. To 
lllustrate, the movement of the target submarine in AHS-l 
during any time interval is represented by equations involv- 
ing the course and speed of the submarine during this in- 
crement of time. Given course, speed, and time of travel, 
the distance treversed and tne direction can be predicted ex=- 
actly. The course itself thouzh, is assumed to be a random 
veriable distributed uniformly on the interval | 0,360) degrees. 

COURSE = (000+RN) 
where RN is a random number drawn from a set of numbers dis- 
tributed uniformly between zero and 360. fTherefore, the di- 
rection traveled during the time interval is a random variable 


anc the todel is stochastic. Once the random number has been 
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draym, COURSE becomes the observed value of 2 random viri2bl. 
and the direction of travel is uniquely determined. 

The word simulation is used interchanseably with model 
by wamy writers in the field. We chceose to @@fine 4 iis 
tinction between the two. As used in this paper, a simula- 
tion includes both a model and its losic. Consider the model 
as a plan which is essentially static until pui into effect. 
The simulation is the dynamic implementation of this plan, and 
includes the program logic and the completed computer instruc- 
tions. When the term simulation is used in this thesis, com-~ 
puter simulation will be implied unless the context makes it 
clear that this is not the case. A computer war same is one 
particular type of simulation. 


&. S. Householder defines the Monte Siarlo method as "..i 








the device of studying an artificial stochastic model of a 
physical or mathematical process..." [ 4] The term Monte Carlo 
ls relatively new but the method is the technique of model 
sampling used by statisticians for at least sixty years. An 
alternative and perhaps simpler statement of the above defini- 
tion is "the use of sampling to estimate the answer to a math- 
ematical problem." B For a comprehensive treatment of 

Monte Carlo techniques and applications the reader is referred 


+0 Symposium on Monte Carlo Methocs, edited by H- aA. Meyer. | 6 | 


The temnct real world and system will be used throughout 





the thesis and in this particular context are considered to be 
synonymous. Both terms are used to describe the particular 
part of the universe which is being studied. The real world, 


or system, could comprise the particles existing in the nuclens 
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of en atem, or it might comsist Orie vortiom of the oceen, a 
submerged submarine and a flisht of soner equipped helicopters. 
In the latter example the weather, sonar operators and pilots, 
or anything else which might interact with the helicopters 

or submarine during a search operation would be part of the 


system. 





CHARTERS IT 


VOB PUTER SIPULATION DSsIcn 


The first principle of computer simulation design is 
that every design situation is uniaue. Each simulation pre- 
sents problems which must be solved on the basis of the ob- 
jyectives of the particular study and the characteristics of 
the real world system under consideration. t is believed, 
however, that certain factors of a general nature must be 
considered in any computer simulation design and that there 
ic 2 locical order in which these factors should ordinarily 
be examined. The purpose of this chapter is to describe a 
general methodologic approach to computer simulation design. 

Before attempting to describe a methodolosy of general 
applicability to computer simulation the question of why such 
a procedure is needed should be answered. War saning has a 
long history of application to military problems and in the 
course of this history many principles of game design have 
evolved and been recorded. Unless there are some basic dif- 
ferences in the objectives as well as techniques associated 
with the desirn process, these principles of classical war 
eamedesicn should be equally applicable to computer simulae 
tion. In an attempt to determine if such differences do exe 
ist, some of the factors contributing to the philosophy of 
both classical war game and computer simulation desirn will 


be considereca. 
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2.1 Effect of the grainime role on classical War @ame (design. 

war games began to achieve the stature of a serious busi- 
ness rather than a form of amusement in the early 19th century. 
2,7] War games from this time on were for the most part exe- 
tensions of the field maneuver. Field maneuvers provided a 
means of practicing for war without the high costs in life and 
resources associated with actual combat. Similarly war games 
offered militarists as well as prospective military officers 
the opportunity to practice the business of war without field= 
ing troops and their associated equipment. 

In both field maneuvers and manual war games, the primary 
purpose is generally training of the participants rather than 
analysis of the effects of factor interactions on the outcome. 
Because people tend to learn more rapidly when placed in fami- 
liar surroundings, the game designer attempts to create a sinm- 
ulated environment which will be familiar to the players. As 
a result, game factors are chosen because of their relation 
to the appearance of reality rather than the effect these fac= 
tors might or might not have on the final outcome of the game. 

If general trends or cause-and-effect relationships be- 
tween numerical or qualitative inputs and game results are 
noted, they are often considered as a secondary return rather 
than a primary fame odjective. Consequently the analysis of 
results, and the question of what specific factors to analyze, 
is gener2lly not considered a part of game design, and is 
often left entirely to the decision maker or his staff. The 
nature of this analysis is also related to the role of war 


gaming as a training device. Throughout the plav of the tame 
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a detailed battle history is usually maintained. This history 
is then studied from the standpoint of the effect of each in- 
dividual decision on the final outcome of the same. If the 
game is repeated, the players and decision rules may chansc. 
In this way the participants are offered the opportunity to 
react to a wider range of stimuli. This is desirable from an 
educational standpoint. However, when war games are used for 
analytic purposes this type of game may test the player rather 
than operational plans or doctrines. 

2.2 Effect of dicital computers o ame desiren, 

The digital computer has contributed more than computa- 
tional speed to simulation. With this means of rapid calcula- 
tion the analyst has gained a degree of control and repeate= 
ability not attainable in manual war gamins. Unlike manual 
raming in which decision rules and players might change be- 
tween plays, these factors do not vary between revlications 
of a computer simulation. Neither will the computer simulation 
results be affected by the participants learnings the same nor 
by some unusually bold stroke of genius exhibited by one of 
the players. Through the elimination of human intervention, 

a major source of uncontrollable experimental error, and the 
gain in repeatability, computer simulation offers the analyst 
an analytic rather than a primarily educational tool. 

The automation of war gaming has also placed certain re- 
straints on the simulation desigmer and these restraints must 
be reflected in any plan for simulation design. The amount 


of information which can be stored in the memory of a computer 





is limited and computer runnins time is expensive, on the 
order of ten dollars a minute for the IBM 7094. FPurtheraore, 
before the simulation can be run on a machine it must be re- 
duced to a set of computer instructions. This list of in- 
structions, or computer program, may represent 2 considerable 
investment in time and money. Large computer simulations such 
as the Air Battle Model described by R. H. Adams and J. Le 
Jenkins | 8 | may be two years or more in the construction. | 9| 
Because of these inherent limitations, it is usually necessary 
to exclude certain factors which do not contribute substan- 
tially to the purpose of the game and may even detract from 
it by masking the same response to factors which do. The 
problem confronting the simulation designer is not just how 
to achieve a high fidelity representation of the real world, 
but how to selectively choose those components of the real 
world which characterize the particular questions to be 
answered. <A major effect of digital computers on game desien 
a¢ thet a need has arisen for criteria which will help 

the designer in deciding what is and is not pertinent to 

the game purpose. 

2.5 Qbjectives of a computer desicn methodology. 

The preceding discussion suggests that computer war game 
ing entails somethinz more than adapting manual sames for vl2y 
on high speed computing machines. It seems resonable then 
that coaputer wame agesign should embody more than a strai¢ht- 
forward extension of manual war game techniques. The fore- 
sOin; also indicates some of the requirements a plan for sin- 


ulation design should satisfy. MThese are; 
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(1) The simulation should be desisned with svecific 
questions in mind. 

(2) The design should capitalize on the advantages of 
repeatability and precise parameter control offered 
by high speed compmters. Tais implies Wreilizow 
some means of reducing tne effect of exverimental 
error. 

(3) mpvhasis should be placed on diminishing the dis- 
advantages of limited space and high cost associated 
with digital computers. 


(4) FPactors should be chosen on the basis of their re- 
lation to simulation objectives. 


(5) Planninz of the analysis and interpretation of re- 


sults should be concurrent with designing the sim- 
ulation. 


2.4 Experimental desis approach to simulation. 

It is considered that the above requirements are best inet 
by the methods of experimental design. This is not surprising 
Since a simulation is basically a sampling exverliment. The 
principle departure of simulation from the agricultural ex- 
periments, wnich save experimental design its impetus, is 
founded in technique and desree of abdstraction rather than 
theory. In simulation the experimentor first builds 2 re- 
presentation or model of tne real world and then performs his 
experimentation on this model. The agricultural worker may 
use real plants grown in real soil as the basis of his experi- 
ment, but not just any plants grown on just any soil. The 
plants are carefully chosen and plots are arransed in 2 se- 
lected patterm in order to control the effect of extraneous 
factors such as soil variability or differences in drainage. 
Similarly, the simulation designer selects factors according 


to whether they contrlbute to the essence of realism in 2 


1 





specific situation. The agricultural exnerimentor eliminates 
weeds which might, by taking water from the plants, result 

in an erroneous estimate of the characteristic yield. A 
problem confronting the simulation desisner is that he must 
first determine which are "vlants" and which "weeds". Once 
this determination has been made he must resist the tempta- 
tion to simulate the "weeds" just because they are part of 
the real world he is dealing witn. 

A further motivation for tne use of the exverimental 
design approach is suggested by I. Fe Warner. 

There have been times when the inexverienced War Gaming 

researcher has set up his own study, planned his own runs 

and after the runs have been made, offered his data for 
someone to interpret. The result is usually an inability 
to place proper interpretation uvon the results, nec= 

essary runs have been omitted, redundant and superfluous 
runs have been conducted and the time and money exvended 

may exceed the value received. [ 10] 

ImpLicit in this~approach is that sé@iection of factors 
to include in the model is closely related to data collection. 
The collection of data is as imvortant a part of statistical 
analysis as data preparation and interpretation. 


For those who may not be familiar with basic experiment- 


al design the followins from Design and Analysis of Industrial 





Experiments, edited by 0. L. Davies, should further clarify 
the relation between experimental design and the discussion 
co Tollow. 


The first step in planning any experiment is to form a 
clear picture of objects of the experiment, the factors 
which may affect the results, and the errors whicn will 
inevitably arise. The chief objects of the experimental 
design are: 


Gig: " arraénge@ the eXp@rimént so that the effhct of 
cnaneing each condition can be readily measured 
and separated from the effects of chanrcin~ the 


De 2 





other con@itions, and from experimental error. 

(Ge) Po obveaiw @ valid estimate of error “appropriate 
for assessing the statistical significance of 
the effects of the factors. 


(141) To enable the effects to be measured with the 
required accuracy. [1] 


The preceding principle of experimental design will play 
an important role in determining the course to be followed in 
designing a simulation. The first step will be to form a 
clear picture of the objects of the experiment and a first 
approximation of the factors affecting the results. This is 
treated in Chapter III, and will take the form of a prelimin- 
ary study of the simulation purpose and the system to be sim- 
ulated. 

Consideration of the sources of error and the first ob= 
jective of experimental design stated above will shape the 
main body of the provosed simulation design method. The nec- 
essity of designing the simulation so that the effects of 
varying parameters can be measured and separated from the ef- 
fects of errors and each other was seen to be a basic departure 
from classical war game requirements. This objective will be 
realized in the computer simulation by designing the model it- 
self as a factorial experiment, and is the subject of Chapter 
IV. The second and third objectives of exverimental design 
are intrinsic to model testing, and will be discussed in con-~ 
jJuaction with implementation of the model in Chapter V. 

Before proceeding to the preliminary study phase a brief 
discussion of factorial experiments is in order. MThis class 


of experimental designs wil constitute the pattern for the 
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simulation model. 
2.5 Factorial experiments. 

Factorial experiments represent a class of experiments 
which are particularly applicable to studying variation de- 
liberately introduced by the experimentor. [11,22] The specific 
techniques associated with constructing these experiments is 
beyond the scope of this thesis and will not be considered. 

The intention here is to consider certain very basic concepts 
of experimental design rather than techniques. The purpose of 
this discussion is to indicate how these concepts can help the 
simulation model designer in selecting inputs which reflect 
the objectives of the system as they pertain to the simulation 
purpose. 

The following terms will be used throughout the remainder 
of the discussion of factorial experiments and simulation 
design. 

Factor. In general a factor is any quantity which cone 
tributes to the determination of the state of the simulated 
system, and can be controlled by the experimentor. Examples 
from simulations of military operations are number of aircraft, 
weather conditions, or electronic counter-measures. 

Factor Levels. The Values assigned to a factor are term- 
ed levels. A level may represent a number such as thirty 
bombers or it may be qualitative. Jamming or not jamming 
would represent two levels of electronic counter-measures, 

a qualitative factor. 
Treatment. A treatment refers to the level of the factors 


Which describe one simulation run. A convenient notation for 
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a treatment is (A, B,0,)- This indicates a simulation run with 


-~ 


factor A and °C at their lower levels and factor B at the second 
level. 

Response. This is synonymous with simulation results as 
we have used this term. A treatment response revresents the 
result of a trial of an experiment corresvonding to that 
treatment. 

Main effect. A main effect for a factor at two levels 
is the average difference in response due to the change in 
levels. 

Interaction. An interaction is said to occur when the 
response to changes in one factor depends on the particular 
level of another factor. An example might be artillery and 
spotters. Increasing the number of firine pieces might not 
increase the damage unless additional spotters were assicned 
to direct their fire. Similarly increasing the number of 
spotters might have no appreciable effect by itself. 

Sample. A sample consists of a series of plays which 
differ only as a result of the series of random numbers used. 
Samples consisting of more than one play are generally used 
in simulaions whose output is either success or failure rather 
than a detailed battle history. When a simulation run con- 
sists of more than one sample the additional samples represent 
replications of the same treatment combination. Replication 
vrovides the means of estimating experimental error. 

The factorial experiment was chosen as the pattern for 


Slmulation design because it is generally the most efficient 
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effort. In ezsente tne design and snalysis of #@ tactorial 
experiment consists of: (1) conducting trials of the experi- 
ment with coubinations of several factors at uore than one 
Level, (2) observing the response to each treatment and 
determining estimates of main effects, factor interactions and 
experimental error, (27) testing hypotheses to determine which 
of the effects are real and which are attributable to e-ver- 
imental error. 

An individual experiment on a completed simulation wav 
be conducted 2c a complete factorial experiment in which ob- 
servations are made with every factor at every one of lts 
Levels. This ugu@lly requires © great many computer runs, 
however, and some form of fractional factorial exnerinent&tion 
might be employed. The eSsential point 1S tnat by buildings 
ai. of the Ta¢etors which are pertinent to the study puppow 
into the model from the very beginning, the cavabilit,y of 
efficiently examining sny combination of these fectors end 
factor levels will exist in the completed simulation. In 
this way additionel assurence that the analysis will reflect 


the study purpose will have been gained. 





CHerien Bl 


PRELIMINARY STUDY 


The design of a computer game, as with essentially all 
operations research studies, must besin with a statement of 
the problem. In most cases the initial statement of the 
problem will come from the person or persons requestine the 
Stucy. The concern here is more with the interpretation of 
the problem by the simulation designer. 

3-1 Formulation of the problem. 

The first step toward formulating the problem is to be- 
come familiar with the various components of that portion 
of the real world from which the problem comes. Components 
in this context may include hardware such as weapon systems, 
Operating doctrine, overational environment and a myriad of 
other items Which make up the particular aspect of the real 
world under consideration. It is suggested that a block flow 
diagran showing the functional position of each component in 
the system be constructed as part of this initial research. 

A flow diagram is also useful for showins interactions among 
components. 

The determination of the scope of the proposed design 
should be concurrent with the preceding step. As the analyst 
familiarizes himself with the system components, the size 
and complexity of the problem area will become more apparent. 
Conversely, unless he has some feel for the rance of activity 
included in the problem it is impossible to know what the cor- 


ponents of the system are. 
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if the Gesian is to be built around varvicilear questtiows 
as related to the oodjectives of tne system, these questions 
must be clear from this point on. It is therefore prudent 
to compare the interpretation of objectives with stated ob- 
jectives. It is possible that on the basis of the above 
initial research, the original problem may be modified or 
expanded, or another research method chosen. 

3.2 Applicability of computer simulation to the vroblem. 

Once the scope of the problem is known and the designer 
hag gained some understanding of the system, it must be de- 
cided whether or not computer simulation techniques fit the 
Situation. This is 2 step which is perhaps too often omitted 
in simulation design. The following statement by C. J. Thomas 
concernins unfortunate practices in milltary gaming is equally 
applicable to computer simulation. 

eeethere is an unfortunately common tendency to think 

of "military gaming" as being not only universally 

applicable but also mechanically applicable. The failure 


to match technique with problem is often expensive and 
almost invariably produces disappointment. [23] 


Consideration of theoretical model. It is generally 
considered that simulation should be used only if it is 
impossible to solve the problem using a simple theoretical 
model. | 14,15,16 | In applying this generality it is important 
that we do not lose sight of the given problem in trying to 
fit the system to an analytic model. Any situation can be 
treated analytically if it is simplified enough. However 
Ooversimplification of the model and the subsequent analysis 
may invalidate the results of the study. Consider the hreli- 


copter search situation treated by AHS-1. The system itself 
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is not overly complicated, and analytic models of a search 
helicopter versus an evading submarine conflict have been 
built. One of the helicopter search plans currently opera= 
tional with the U. S. Navy is based on a game theory solution 
of this conflict situation. The model used in this solution 
includes the assumption that submarine maximum speed is known, 
and that helicopter sonar ranges and dip cycle time are known, 
and the same for all helicopters. The simulation study, AHS-1, 
was motivated by the fact that these factors among others are 
subject to uncertainty. The author wanted to know the rela- 
tive effectiveness of certain search plans, including the so 
called optimum plan found analytically, in the face of this 
uncertainty. 

Expected number of plays. A second consideration in de- 
Cciding when to resort to computer simulation, is the number 
of expected repetitions of the experiment. The validity of 
the Monte Carlo technique is dependent on repeated trials. It 
is not unheard of to make ten or more consecutive sevens with 
two fair dice. Inferring that this is routine on the strength 
of observing the phenomenon once, however, could be financially 
disastrous. Similarly, confidence in simulation results is 
related directly to the number of replications. 

The number of input parameters to be examined and the 
range of interesting values of these parameters also affects 
the number of required repetitions. The effect on experimental 
results of varying parameter values is usually of more interest 
than the particular outcome corresponding to any one set of 


input values. In addition certain parameters may exhibit 
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interaction. That is the effect of changing the input value 
of two or more factors simultaneously is considerably dif- 
ferent from that of varying these same parameters individ- 
ually. 

It is apparent that if several parameters are to be ex= 
amined at more than one level, each run consisting of many 
replications, the number of simulation runs can become very 
large. Even with a game of relatively limited scope such as 
AHS-1, several months might be required to examine all input 
parameters at only a few levels by manual methods. On the 
other hand, if only a few plays of a game are contemplated, 
computer simulation is probably not the most efficient method. 
A completed computer program represents a sizeable investment 
and the cost in time and money of running this program on a 
computer one time represents a small fraction of this invest- 
ment. 

Player training. The potential value of training derived 
from player participation in the game must also be considered. 
If the primary purpose is clearly player training then a 
manual game is indicated. Whether the purpose is educational 
or analytic is not always tacitly obvious, however, and a 
decision to satisfy one or the other requirements will have 
to be made. 

Assuming that computer simulation is theoretically ap- 
plicable the analyst must next ask if it is physically prac- 
ticable. A study which is worth doing is worth doing correct- 


ly. If the limitations of time and available phisical faci- 
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lities preclude an adequate treatment of the problem by sim- 
ulation methods the person requesting the study must be ap- 
prised of this. 

3.3 Research into existing simulations. 

Having decided to use simulation methods, the designer 
might next ask if there is an existing model which is suited 
to this particular situation. However, in many respects try- 
ing to fit the problem to an existing model is contradictory 
to the stated objective of tailoring the simulation to the 
particular study purpose. It would be preferable to re- 
state the question more generally as; what is the existing 
level of simulation knowledge which may be related to this 
system? Although it might be unwise to modify the particular 
situation to fit an existing model, it is equally foolish not 
to profit from the effort of others whenever possible. 

fo illustrate this point we refer to the scheme by which 
sonar detection ranges are computed in AHS-l. The particular 
objective required some more realistic means of determining 
detection ranges than the usual procedure of using a determ= 
inistic assured sonar range. The scheme selected was design- 
ed by R. L. Klinkner as a general method, applicable to most 
existing Navy sonars. {17| Furthermore an alcorithm for con- 
Duting tne needed sonar ranges by this method had been writ- 
ten and coded in the computer language previously selected 
for AHS-l1. The availability of this computer coding made it 
possible to design a simulation which fit the requirements of 
the problem within the time allotted. This misht well have 


been impossible otherwisc. 
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CHAPTER IV 


DESIGH "OP THE MODEL 


Tges model of inwverest im thise chapter is a@ stochastic 
mathematical model. Recall that a mathematical model was 
defined as a mathematicai representation of some aspect of 
the real world. This model can usually be subdivided into 
models of each of the individual components or activities 
Which characterize the cystem of interest. It is partly 
because of this subdivision that the construction of a block 
flow diagram was susgested as vart of the formulation of the 
problem. Such a flow chart is often a useful aid to recos- 
Hition of bogical functional subdivistone. 

4.1 Functional components of the model. 

There are several reasons for subdividing the model in- 
to functional comvonents, the degree to whicn this is done 
beinz closely related to the scove of the system under study. 
4ny system, even one of limited complexity such as the helicon- 
ter-submarine conflict of our examole, is much simpler to 
visualize if it is broken down into functional components or 
activities. 

A second arzunent for subdividing tne model alons func- 
tional lines is that for more complex problems several versons 
may be simultaneously engaged in modeling the system. Con- 
sider a simulation of a large area anti-submarine operation. 
Helicovter operations might constitute one logical subset 


to be allotted to one individual. One person could develop 
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this part of the simulation, code it and complete the vrocran 
check-out independently of the main study group. The example, 
AHS-1, could with only minor modification represent such a 
subsection of the larger simulation. 

Dividing the model into such subsets also facilitates 
parameter validation and program check-out. This applies 
equally whether the simulation is being built by ome person 
or several. 

The preceding model subdivisions are primarily for mech- 
anical or operational convenience of the designer. They are 
aids to organizing the existing knowledge about the system so 
that the model will become physically manageable. Models may 
also be subdivided into conceptual, or structural, components. 
4.2 Structural model components. 

The structural, or conceptual, components of a model are: 

(1) Assumptions 

(2) Input parameters 

(3) Input Yerbebles 

(4) Algorithms 

(5) Measures of effectiveness 

It will be simpler to define the above terms and explain 
why they are considered to be useful to the model designer if 
the following two definitions are asreed upon. These defini- 
tions are from Van Nostrand's Scientific Encyclopedia. [ 18] 

Parameter. An arbitrary constant, as distin&uished fron 

a fixed or absolute constant. Any desired numerical value 

can be siven to a parameter. 

Variable. aA quantity as distinguished from a constant 


to which any number of values in a siven set inay be assion- 
ed. If thé sét comprises 4 dommin(a,b), then the variable 








is only defined over tnis interv2l and all values out- 
side of the interval are ignored. 


These definitions will apply whenever parameter or variable 
is used in the succeeding discussion. 

Assumptions. Any conclusions arrived at on the basis 
of simulation results must be considered invalid unless they 
were made with a complete knowledge of the major model assump- 
tions. The model assumptions reflect more than any other 
factor the designers understanding of the system to be modeled 
and his interpretation of the objective of the simulation. 

The assumptions reflect what is known precisely about the real 
world as well as what is uncertain and the degree of uncertain- 
ty. To illustrate, the following are model assumptions from 
AHS-1. 

The ocean is represented as a rectangular coordinate 
system in the simulation. This is based on the assumption 
that the area encompassing a helicovter search is small enough 
that the earth's curvature need not be considered. Helicopters, 
which are characterized as points from which associated sonar 
ranges are measured, are assumed to move about this square 
e@rid in discrete jumps. These assumptions reflect certain 
well defined facts about helicopter sonar search operations. 

It is further assumed that the above sonar ranges vary 
among helicopters according to some probability distribution, 
ed that this distribution is a function of sonar figure of 
merit. It is known that sound propagation conditions, at 
any one point in the ocean, are a function of a time correlated 


random variable. | 17,19 | This time dependence of sonar 
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is only defined over tnis intesvil and all values oute- 
side of the interval are ignored. 


These definitions will apply whenever parameter or variable 
is used in the succeeding discussion. 

Assumptions. Any conclusions arrived at on the basis 
of simulation results must be considered invalid unless they 
were made with a complete knowledge of the major model assump- 
tions. The model assumptions reflect more than any other 
factor the designers understanding of the system to be modeled 
and his interpretation of the objective of the simulation. 

The assumptions reflect what is known precisely about the real 
world as well as what is uncertain and the desreec of uncertain- 
ty- To illustrate, the following are model assumptions from 
AHS-1. 

The ocean is represented as a rectangular coordinate 
system in the simulation. This is based on the assumption 
that the area encompassing a helicovter search is small enough 
that the earth's curvature need not be considered. Helicopters 
which are characterized as points from which associated sonar 
ranges are measured, are assumed to move about this square 
grid in discrete jumps. These assumptions reflect certain 
Well defined facts about helicopter sonar search operations. 

It is further assumed that the above sonar ranges vary 
among helicopters according to some probability distribution, 
and that this distribution is a function of sonar figure of 
merit. It is known that sound propagation conditions, at 
any one point in the ocean, are a function of a time correlated 


random variable. 117,19 | This time dependence of sonar 
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ranges has not been included in the model for two reasons. 
First, a helicopter search occupies a relatively short time 
span and the variation due to the aoove random variable is 
expected to be small in the time interval of interest. 
secondly, sonar ranges can be expected to decrease or in- 
crease with the same relative frequency during a search. If 
for some reason sOnar ranges were biased so that either a 
decrease or an increase with time could be expected this 
assumption might not be valid. 

The assumptions concerning sonar range all have a bear- 
ing on uncertainties in nature and the effect of those une- 
certainties on the simulation purpose. All of the preced= 
ing assumptions reflect the model builder's understanding 
of basic functional relationships of the system components 
as they are related to the simulation purpose. 

Input parameters. These are constants which may be 
assigned any arbitrary value by the simulation user. Our 
interest in input parameters is generally in the way changes 
in parameter values affect the outcome of the operation sin-~ 
ulated. The factor which characterizes these quantities as 
parameters rather than variables is the derree of human con- 
trol in the real world situation being simulated. For example, 
the varticular search plan designated by a helicopter flisht 
leader is entirely within the control of the flight leader 
in question. In assigning values to this particular input 
parameter the analyst is concermed with finding criteria by 
which a certain plan should be selected, given a specific 


state of nature. Note that input parameters need not be 
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nmamericS) quantities nor Will they necessarily of Lhe actual 
inputs to the ororram. Search plan is 2h inpwy peraneTer 
and it® vVelue mitht be Bhan Lisa or Snirel Seewen three. 
This quantity will givé rise to invut auentitieés Which do 


neve numerical values, however. In this ec ample bdbeartines 
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and distances of tne serrch legs are such numerical valued 
quantities. 


Input variables. Tnese ar 


we cannot expect to have control in actual vractice. For 


{D 


input nvantities over which 


example submarine speed, sonar conditions or tine late at 2 
Summerine daium. Inpwt variables ere chamacterigtic of "ie 
system and the ranse of values assigned to these auantities 
is 2ssenilally determined by nature. 

As noted above for input parameters, input variables as 
initially defined ty the model deeisner may not be the actual 
Drozratza inputs: Ih most cases tnese Will be determined im 
conjunction with the detailed investiration and data collec- 
VLOm intrinsic to tite aciwal model comsirnuetlon. 

To tie the above definitions toctether, consider the func- 
clomel relation from MAis-1: 


a (ST 52 Bers) 


a’ N 


Output (Number of de‘Vections during a run) 


wnere 


O 
5. = Submarine Speed 

< Inout variables 
tT, = Time late of helicopter at Gatumr 
Pp. = Searen Plan ¢€ 

A Input paraneter.: 
H,. = Sunder of helicopters 

Is 





F is the rule, determined bv the model assunptions and 
simulation logic, Which determines the dependence of Ry 
tne input variables and parameters. Submarine speed and tine 
late are independent variables which are allowed to vary over 
selected portions of their assumed ranges. The input pvara- 
meters may be assigned any arbitrary set of values. Whereas 
the range of input variables depends on what the analyst has 
assumed to be the possible states of nature, input parameters 
reflect particular questions about the system. The optinunm 
values of these parameters corresponding to specific states 
of nature, i.e., specific values or ranges of values of the 
input variables, revresent pnart of what we hope to learn 
from the simulation. 

Input parameters and variables may change roles accord- 
im@ to the specific circumstances of application. For 
instance, it may be possible to control time late at datum to 
some extent but at the expense of the number of helicopters 
Which can be sent to the datum. It might happen tnat two 
helicopters are airborne when a submarine datum is established. 
The decision minust be made to launch more helicopters and 
accept a sreater time late, or proceed immediately to datun 
With the available aircraft. Here time late assumes the role 
Of @ parameter and number of helicopters is essentially a 
variable. 

Aleorjthm. "4n algorithm is stated rule or procedure 
that may be followed to arrives at some specified coal. "[20] in 


algorithm is characterized by the fact that 1% tay be 


—_- aca 
— ——_- 2 >a 
_ °° === om eo 
a &2eMese ia eae” ae 
= iad’ . = see Cop @ © 


—. — ‘. o- * Gummi Aaeae 
= —] EE ee es —& 


® ea «* zs «& qe &> —eEeE~.U—C<iC SoS] 


—_——> 





mechanically applied. ‘he procedure for seneratin® rendcom 
numbers employed in aHS-1 is an algoritnr. 

The model designer must consider every vossible question 
Which might arise and preprogram an answer to each question. 

It 1s not possible to defer a decision until such time durines 
the play of the game that the need for making the decision 

may arise. Therefore, in order to build a model which can be 
incorporated in a computer program it must be possible to 
reduce the essence of the real world situation to a collection 
of algorithms. By means of this set of algorithms the program 
operates on input variables and parameters in a way determined 
by the assunptions. 

Measures of effectiveness. Saaty defines measure of ef- 
fectiveness as "a criterion which measures the extent of suc- 
cess of a solution, as related to the objectives, when 2p- 
plied to a problem arising in an operation." [22 | The prin- 
ciple measure of effectiveness incorporated into AHS-1 is 
number of submarine detections. In this example the objective, 
as used in the above definition, is to evaluate the comnarative 
effectiveness of various search plans in similar tactical 
and physical environmnents. One of the problems in this ex- 
ample is that the number of available helicopters is generally 
limited. 

4,3 


eterminine measures ectiveness,. 





Military simulation studies fall into two broad cate- 
s0riés according to purpose; testing vilans and evaluation of 
weapon systens., Regerdlees of where the simulation fite in 


hie @wencral clateifive@tion, however, if decisions are to be 
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influenced by tne simulation results some measure of the suc- 
cess or failure of the operation in a given instance is re- 
quired. Measures of effectiveness may be comparatively 
obvious, or quite subtle. In the case of AHS-l1, the number 
of submarine detections is a straightforward measure of the 
comparative effectiveness of two search plans under similar 
circumstances. 

Consistent well defined measures of effectiveness are 
essential to good simulation design. To a considerable extent 
the measure of whether an included parameter or variable is 
essential to the simulation purpose is the sensitivity of the 
measures of effectiveness to changes in this input. vonsequent- 
ly, the measures of effectiveness will represent imvortant 
criteria by which the analyst must decide whether a factor 
1s to be included or omitted from the simulation. 

Inconsistent measures of effectiveness indicate incon- 
patability of the simulation objectives as well as the lack 
of a clear interpretation of the objectives of the system 
being simulated. In essence the measures of effectiveness are 
the statements of the objectives of the system beins studied. 
The selected measures of effectiveness are indicative of the 
designers understanding of the system objectives just as 
model assumptions reflect his understanding of basic func- 
tional asnects of the systen. 

Output data. As far as 1s possible the output of the 
Simulation should be determined in conjunction with selecting 
measures of effectiveness. All of the measures of effective- 


ness will be output auantities and these may also sive rise 
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to associated numbers which will be of interest to the 
analyst. In addition, certain control variables should be 
included in the output. These controls should be quantities 
whose values or distributions are well enough known that they 
can be used as rough checks on the logic and numerical inputs 
to the simulation. 

4.4 Selectinze inputs 

Everything that has been accomplished up to this time 
has been preliminary to the actual design of the model. During 
the formulation of the problem the scope of the simulation was 
determined. The scope in turn determined the accentable 
level of aggregation consistent with the stated problem. Hav- 
ing chosen measures of effectiveness by which the simulation 
results, and to some extent the simulation itself, will be 
evaluated; the designer is prepared to select the inputs and 
incorporate these factors along with certain assuwnptions and 
rules into a representation of the systen. 

Inputs were classified as either variables or parameters 
on the basis of the extent to which these quantities can be 
controlled in the real world. This classification will be 
useful in determining the ranges of input variables and para- 
meters and the distinction should be kept in mind durines the 
following discussion. In determining which inputs to include 
in the model, however, it is more suitable to our purpose to 
combine inputs and variables under the single title of factors. 

D. R. Cox lists four items to be considered in desisnine 


ae Factories. experiment. | 12| These are, the factors to 
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include in the design, the levels of each factor, the number 
of observations to be made, and measures to reduce the effect 
of uncontrolled variation. The last two of these are more 
closely related to simulation testing and analysis of the 
Simulation results. The factors to be included and the levels 
of each factor determine the form of the model itself and will 
be considered at this time. 

It should be noted that at this point the simulation 
designer is concerned with a general experimental design. One 
particular experiment might involve only two or three factors, 
each at one or many levels. However, the simulation objece- 
tive often involves several questions about the system, and 
the analyst must design the simulation to encompass a class 
of experiments. Consequently, all of the factors necessary to 
conduct any experiment of this class must be built into the 
model. 

Choosing factors. Whether a particular factor should be 
included in the simulation is primarily a matter of exverience 
and judgment. Once a working version of the simulation has 
been completed, statistical methods can be employed to test 
the sensitivity of the measures of effectiveness to the select-~ 
ed inputs. In addition a certain amount of parameter testing 
can be accomplished by hand or desk calculations as the model 
is being built. For the most part, however, factors must be 
selected on the basis of experience, that of both designer 
and user, conditioned by what was learned durins the preline- 


inary investigation of the problem. 


Bak 


> — i> —ia, = 
eso 7 == ho c 
-_¢ &t¢ @ lee 


a. 


es iy mn ene —— ee 
= —_— = mii ie —z 
SS EE OO Om 
~ at ati—uii—, da =_ => 
— ————, © 7 
—eEeEeEeEeE 


= 





An example from AHS-1 may lend more meanins to this 
assertion. Probability of gaining a non-submarine contact 
and maximum time to evaluate such contacts are included as 
input parameters. Actually these represent anv delay exper- 
jenced by a helicopter while conducting a sonar dip; non-sub- 
marine contacts are one of many reasons for such delays. 
That delays occur is well known to those who have pariicivated 
in helicopter search overations. However, the existence or 
significance of variation in dip times would not necessarily 
be apparent to an analyst who had not varticipated in anti- 
submarine overations, solely from a review of the literature. 

This introduces another question to be considered oy the 
model desigmer. To what extent may too close association with 
the day to day operation of the system limit the persrective 
of the simulation customer? For example, a catapult officer 
may be a valuable source of information concerning aircraft 
carrier air operations. However, this officer is involved 
with many factors which are essential to air operations but 
may not be pertinent to the purvose or level of acereration 
of a simulation of air operations. Gaining the maximum in- 
Pommetion: frem this officer, Without clutterimg the Simuletion 
with unnecessary detail also based on the man's exverience 
is again a matter of judgment. 

The designer must also appreciate the role of the decision 
maker when considering inouts. It is the user wno rust Lake 
the final responsibility for decisions based on ite sinulea- 


tion results, and hence these decisions will be influenced to 
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a larse extent by his confidence in the simulation. In the 
above exampvle if exclusion of catavult pressure from the sim- 
ulation will result in a complete lack of faith in the simu- 
lation by the decision maker, then catapult pressure may be 

a valid input even if the simulation results are insensitive 
to this input. It may be that only by including this factor 
and conducting a more detailed analysis of the sensitivity of 
response to it can the question be resolved to either the 
designer's or user's satisfaction. 

Factor levels. In selecting factor levels, just @s with 
the factors themselves, the designer will usually want to in- 
Clude a certain amount of generality in the model. The factor 
levels chosen at this time will in some cases be the particu- 
lar levels used for the analysis later on. In most cases, 
however, the designer will be interested in a ranse from which 
specific levels may be chosen at the appropriate time. In this 
vnase of simulation design, planning the experiment around 
specific questions can be varticularly productive. If only 
certain selected ranges of factor levels have a bearing on 
the study objectives, time spent determining the entire ranze 
of possible values might be essentially wasted. Restricting 
the levels to interesting ranges can also result in many ef- 
ficiencies in the subsequent computer program. das the range 
of factor levels is \Increased so is the computer storage 
Space necessary to contain the associated numerical values. 

Determining factor levels will also sive rise to the 


need for determining certain distributions of associated 
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mput variables. It is in this stev™that the data collection 
Pole of Statistical analyset@ is of obrticuler interest’ to 
the simulation desirner. 

Input parameters are directly related to the particular 
questions vosed by the person or activity requesting the sim- 
ulation study. Consequently the customer is often the best 
source of the range of parameter values. Input variables, 
however, are more closely associated with the possible states 
of nature. The designer is free to include or reject a speci- 
fic variable according to its pertinence to the simulation 
purpose. Once he has included a particular varlable in the 
model, however, the analyst must determine tne rance of values 
of the variable according to what actually exists in the real 
world. This determination might involve a comvaratively 
Simple library search or an extensive analysis. Determination 
of factor levels may at times give rise to analytic or simu- 
lation studies which are more extensive than the original 
simulation. 

It might happen that the reauired information cannovt be 
obtained in the time allowed or at a cost consistent with the 
designer's operating budget. Pernavs the information is not 
available at any cost, or only a very restricted sezment of 
the true range of the variable can be determined with any de~ 
eree of certainty. The designer may leave this variable out 
or he may resort to an educated guess. ‘What is done Js a 
matter of judgment and depends on the particular situation. 


It is essential, however, that the simulation reflect tne 
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nature and degree of uncertainty relative to this input. 

As with any analysis the results nave validity only with 
respect to the assumptions and input data from which they come. 
4.5 ‘Combinine factors and assumptions in the model. 

As an illustration of how the assumptions, factors and 
factor levels, and algorithms specify the model, we will con- 
sider sonar detection ranges in AHS-1. From the preliminary 
study it was determined that the simulation purpose dictated 
that sonar detection range was an essential element of the 
helicopter model. This quantity 1s essentially dependent on 
nature and was therefore defined as an input variable. 

Further investigation revealed that detection ranges de- 
pend on certain variables such as figure of merit, water con- 
ditions and sea state. These were tentatively assigned as 
factors and the process of determining, general factor levels 
was commenced. As noted above, assigning factor levels is 
closely related to determining the range of input variables, 
and data collection. 

The first step in assigning general factor levels was 
to determine the functional relationship between the orisinal 
input variables, sonar detection range, and the tentatively 
selected factors. Investigation in this area revealed the 
stochastic nature of figure of merit and resulted in a 
modification of the selected inputs. That is the distribution 
parameters, average figure of merit and variance, replaced 
figure of merit as inputs to the simulation. This of course 


involved certain assumptions which have already been mentioned. 
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On the basis of the above research and assunptions the 
valid ranse of factor levels was determinec. It was then 
only left to discover algorithms, by which the specific 
values of the original input variable could be computed, and 
combine these with other algorithms which would move the 
helicopter model over the plaving area. t was of course 
necessary that these algorithms not only be mathematically 
correct but that they be avplicable over the ranse of both 
factors and the original innut variable. The end product 
of this process was a point which could be moved about a 
square grid, and an associated detection range. This is one 


model of a helicopter. 
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5.1 <Addine losic to the model. 

The model is the heart of the simulation, but to imple- 
ment the model it is necessary to move the simulated system 
components throucsh time as well as space. The simulation logic 
is the means by which the model, a static representation of 
the system components, is transformed into a dynamic time 
ordered simulation of the system. In practice the design 
of model and logic may be concurrent, but conceptually there 
exists a distinct interface between then. 

To illustrate this, we refer to the helicopters in AHS-l. 
As stated earlier, a helicopter is represented by a point on 
a rectangular crid torether with an associated sonar detec- 
tion range. These, combined with a suitable algorithm for 
chansing the crid position of this point, complete the model 
of a helicopter. The implementation of this model in a sin- 
ulation of a helicopter search operation requires some pro- 
cedure for causing certain actions to take place at specified 
times. Each helicopter must move to successive dip stations 
at prescribed times. Having arrived at a dip station, a 
search for the submarine and a determination of the outcome 
of that search must be mades To complicate this picture, 
other helicopters are simultaneously acting and interacting 


with the target submarine. 
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The determination of which actions ars to take place 
and their relative order of oceurrence Could be accomplished 
manually, usine the 2bove model. Similarly, loriecal struc- 
tures other than the one used in AHS=-l1 could be employed to 
implement this same model in a computer simulation. 

The two most frequently encountered logical structures 
are the event store, employed in AHS=-1, and the time-step. Both 
of these structures are essentially techniques for apvroxi- 
mating time, which is continuous in the real world, by » 
series of discrete steps. In an actual operation certain 
events take place either simultaneously or in an overlapping 
time space, but a computer can only perform one overation 
at a time. Consequently when an action occurs in the sin- 
ulation, time must be arrested for all units not participatine 
in that action. ‘Some adjustment must then be made for initer- 
actine events which occur simultaneously. This necessil- 
tates some predetcrmined lorcical order in which actions and 
participating units are to be considered during the play. 

Time-step. In a time-step simulation the play takes 
pi@ee €@8 & series of discrete time increments. The length of 
each interval is called thetime step and may be thought of @2s 
representing the lensth of one fame move. 

During each time step the positions and capabilities 
OF BAA wnitts ere Boted and a list of posSible acthons is 
consulted. For each entry in the list of possible havpenin«s, 
all units which could be affected by this action are then con- 


ef#dered in turn. With respect to each of these units, 
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calculations 2re performed to see if the action in auestion 
did in fact occur and to wnat extent tne wmit Was affectec. 
after every action in ine list has been considered, the 
positions and cavabllities of all particivants are reassessed 
and the play moves into the next time interval. The end of a 
play may be sisnalled by the occurrence of some predetermined 
action or at the completion of a snecified number of time 
steps. 

It if characteristic of this method that @11]1 ewents 
which occur within the same time step are considered to have 
taken place either simultaneously or in the order in which 
they appear in the list of possible actions. For example, 
two possible events might be weapon detonations and missile 
launches. If missile launches are considered ahead of de- 
tonations, then it would be possible for a missile to launch 
when in fact the missile site had been destroyed by a weapon 
detonation prior to scheduled launch time. The nrobability 
of this happening can be decreased by decreasing the lenrth 
of the time step, but at the cost of increased computer run- 
ning time. 

Event Store. aAn event store simulation is played as a 
sequence of events rather than discrete time intervals. In 
contrast with the time-step method, these events are not con- 
sidered until they are scheduled to occur. The scheduling 
of future events is accomplished by entering elements, com- 
monly referred to as event Words, in a table called the event 


store. “ach event word consists of the tyne of event, the 





time the event is to take vlace and any other information 
which is required for the execution of tne event. In the 
sample vrogram, AHS-1, the identity of the unit involved in 
the event is vart of this word. In many simulations, events 
once entered in the list may be cancelled as a result of 
subsequent action, necessitating the inclusion of a validation 
signal in each event word. Events which have been invalidated 
may then be discarded when they reach the top of the event 
store. 

It is characteristic of the event store logic that 
sections of the program representing events are essentially 
independent of each other. Interactions among these separate 
components are caused to take place by two events usually 
referred to as Take Next Event (TNE) and Store New Event (SNE). 
Store New Event enters event words representing future events, 
in the event store in chronological order, and Take Nex: 

Svent causes events to be executed at the appropriate time. 

At the beginning of a play certain events which can be 
predicted from the program logic and the input data are enter- 
ed in the event store. Control of the program then passes to 
TNE which notes the event type of the earliest valid event 
word and calls the appropriate event subroutine into action. 
This cvent subroutine carries out its preassigned function 
according to the information contained in the event word and 
the input data. As a result of this action other events may 
arise, in which case SNE is called upon to store the necessary 
event word. At the completion of each event, control is re- 


tumned to TNE and the cycle repeated. The play terminatcs 
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when some particular action occurs, such as a submarine de- 
tection in AHS-1, or when there are no more events in the 
event store to be vrocessed. 


Time-stev versus Event store. Efficiency of both pro- 





gramming, and computer operation are usually prime considera- 
tions in choosing the logic for a simulation. The event store 
method is somewhat more compatible with the functional sub- 
Givisions of the model previously mentioned. Hach event 
subroutine usually corresponds to some particular function 

of a system component. On the other hand, the vrogrammer may 
have difficulty understandings the chronological sequence of 
events at first. [22] 

Since both structures are methods of approximating real 
continuous time by discrete steps, the number and hence the 
length of each step is closely correlated with machine run- 
ning time. With the time-step procedure one time step may 
be considered as a single approximating step. Every possible 
action and every unit is considered once each time increment, 
60 as the number of time steps is increased, machine running 
time becomes greater. Decreasing the number of operations 
by lenethening the time step results in greater time agrre- 
mation which is usually undesirable. In general the time-step 
method results in greater computer running time when the exact 
order of occurrence of closely spaced events is a critical 
consideration. 

In the event store method each event may be considered 


as one approximation. With this structure, actions besin in 
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the proper time sequence but in the machine only one event 
can occur at any one time. Consequently interactions among 
events which occur in overlapping time intervals may be lost. 
An example of this is the relation between submarine maneuvers 
and submarine detection events in AHS-l. aA detection event 
consists of one or more sonar searches. If the submarine is 
due to maneuver while a detection event 1s in prosress, then 
the correct submarine position will not be available from the 
time the maneuver is scheduled until the completion of the 
event being processed. The effect of this could be minimized 
by definins each individual search 4s one event. This would 
essentially increase the number of events and would result in 
increased machine running time. The procedure that is used 
in this case is to interrunt any detection event between 
searches if a target maneuver is to occur. This is determ- 
ined by asking if the submarine maneuvered at the end of 
each search. The submarine is then maneuvered and the detec- 
tlon event resumed at the time it was interrupted. 

The above is an example of superimposing a time-step on 
One portion of the event store logic. Similarly an event 
store may be incorporated into a basic time-step structure. 
This technique is employed in STAGE a global atomic exchange 
model built by Technical Operations Incorporated under contract 
to the Us, S. Gir Porce. The SRAGE Simus@tor, which is pleyed 
in fifteen minute time intervals, uses the output of a pre- 
processor for its input. Part of the STAG®= preprcecessor out= 


put is a list of events which have positive pronability of 
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occuring during each time interval. This list 1s read from 
magnetic tape into the main simulation program as it is need- 
ed and essentially represents an external event store. fThe 
simulation itself then "Steps through each time period, con- 
sidering the various events scheduled to occur in each time 
period and determining which events actually ocour."| 23] 

Flow charts of the logic. Block flow diagrams can be 
very useful in organizing the logic as they were in delineating 
the functional subdivisions of the model. Flow charts also 
provide a convenient vehicle for communication between simn- 
ulation designer and programmer. 

Flow charts at this stage of the design will usually be 
constructed on at least two levels of detail. General flow 
charts consisting of boxes enclosing plain English statements 
are one means of describing the model and the logical way 
in which model components interact. General flow charts, 
along with the documentation of the structural model compon- 
ents may completely specify the model and logic. However, 
more detailed flow diagrams must usually be constructed be- 
fore a professional programmer can reduce the model and logic 
to a useable computer progran. 

5.2 Computer languages 

Having reduced the model and logic to detailed flow 
charts the designer is ready to turn his creation over to a 
programmer, or as is sometimes the case to write the computer 
program himself. A computer program is essentially a list of 
instructions written in a format, or language, which is use- 


able by the computer. The program is the vehicle by which a 
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complicated mathematical model is reduced to the very limited 
set of operations which can be performed by a computer. 

aA complete hierarchy of computer languages is available 
and the selection of one of these depends on several factors. 
Some of these considerations are: (1) simulation size, (2) 
expected number of runs of the completed simulation,(3) pro- 
gramming facilities, (4) available time for programming and 
program testine, and (5) the available computer facilities. 

Computer lansuases may be classified accordins to re- 
lative programming ease. In general the simpler the progranm- 
mins the more inefficient the program is in its utilization 
of operating time and memory space. One possible classifica- 


tion is as follows: 
Ll. Binary code 


2. Assembly lansuare 
3- Compiler language 
4. Simulation oriented language 
Binary code. This is the only language the computer 
actually "understands" and all other language systems must 
ultimately reduce the programmers instructions to the binary 
number system. Writing binary instructions is extremely 
time consuming since complicated operations must be reduced 
by the programmer to a sequence of basic operations which the 
computer itself can accomplish. These basic operations con- 
Sist of storing and transferring information, addition, sub- 
traction, multiplication and division. 
Assembly or syinbolic lanzuazes. These languages, such 
as IBM FaP and the symbolic language used with the CIC 1604 


employ wnat is called an assembly program. The programmer 
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writes assembly language instructions which are then trans- 
lated into machine instructions by the assembly program. The 
programmer might write a statement such as FAD (R), which is 
symbolic for Floating Point Add (Location 8). This state- 
ment instructs the computer to note the contents of a loca- 
tion previously designated as R and add it to what ever has 
been placed in a section of the machine called the A register. 
This sum will then be left in the A register and may be stored 
in another location by additional instructions. 

This is a much simpler instruction to write than the 
corresponding binary configuration. In addition, the nunber 
of the location in which R is stored need not be Known to the 
programmer, as the assembler will keep track of R once it has 
been defined. When binary machine language is used the program- 
mer must keep track of every location and its contents by 
number. 

Compiler lanzuazes. The use of languages such as FORTRAN 
or ALGOL further simplify the programmer's task. In general, 
however, the step from assembly to compiler language involves 
a considerable loss in storage space and flexibility, and an 
increase in operating time. The FORTRAN program is written 
as a sequence of arithmetic statements such as R = Read. This 
means take the contents of the location assisned to R, add 
this to the location assigned to A and place the answer back 
in R. Here the format is again floating point but this is 
Signaled to the compiler by the particular letters used. 

The FORTRAN program, consisting of statements such as 


the one above and many other instructions which are a freat 
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deal more complicated, must first be compiled. This consists 
of a compiler program translating the FORTRAN statements in- 
to an assembly language program. This assembly language 
program is then converted to machine instructions. This gen- 
erally involves inefficiencies both in the utilization of 
time and storage space. The compiler, lacking the ability 

to reason, is not as sophisticated as a human programmer and 
tends to waste instructions to some extent. 

Compiler languages do have some very important advantages 
in addition to ease of programming. One of these is the large 
library of function sub-routines and programs which is avail- 
able. For instance the simple statement R = SINF(THETA) can 
be used to compute the sine of an angle THETA and store this 
value in the location R.- To accomplish this using an assembly 
language, the programmer would have to write the instructions 
for computing the Taylor's expansion of sine THETA, add the 
terms out to some desired accuracy and place the sum in lo- 
cation R. 

Simulation oriented languages. Languages designed pri- 
marily for simulation, such as SIMSCRIPT or MILITRAN are cur- 
rently being developed. | 24,25) These languages, or program- 
ming systems, are designed to further simplify programming. 

In general, the program is written as a series of functions 
comparable to the event subroutines in AHS-l, rather than 
arithmetic statements. The source program consists of phrases 
Similar to those commonly used in general flow diagrams. 

These statements are then translated into subroutines written 


in a core language such as FORTRAN. 
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Choosing a computer lansuage. Probably the first con- 


sideration in selecting a language is size of the resulting 
program. If the problem will not fit into an available 
machine, running time and programming considerations become 
academic. If the simulation is large, however, the program 
size should be considered in conjunction with available data 
processing facilities. In many cases the effective machine 
size can be expanded far above actual size through the use of 
satellite equipment such as magnetic tape drives. Parts of 
the program and input data may be stored on magnetic tape 
external to the machine. The taped information is then 

read into the computer only when it is needed. 

If space reauirements indicate that a compiler language 
is feasible then the balance between machine operating cost 
and programming expense is very important. If a simulation 
is to be run only a few times, a compiler language is prob= 
ably indicated because of the lower programming cost. Of 
course if programming in assembly or machine language would 
delay the completion until after the information was no long- 
er of interest, a compiler program might be the only altern= 
ative no matter how many simulation runs were anticipated. 

The availability of existing subroutines which may be 
useable in the program should also be given careful consider= 
ation. In selecting a language for AHS-l, the programmer was 
obliged to choose between two compiler languages, FORTRAN 
and NELIAC. This restriction arose because the program was 
intended as an example and these two languages are the more 


commonly used within the United States Navy. The final 
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decision to use FORTRAN was greatly influenced by the avail= 
ability of a sonar range prediction program written in FORTRAN. 

One of the most important drawbacks to FORTRAN for sime= 
ulation studies is the inability to pack more than one quane 
tity into each word. For example the event store table used 
in AHS=-1 consists of three arrays. In assembly language these 
could be combined into one table by allotting only part of a 
word to each of the three items of information. It is possible 
to write packing and unpacking routines in assembly language 
and combine them with a FORTRAN program. Such routines take 
time to program and it might be almost as convenient to write 
the entire code in an assembly language. 

The various dialects of ALGOL, such as NELIAC or JOVIAL, 
do retain the packing feature of an assembly language, but 
at the cost of an increase in computer operating time. 

5-5 Program check-out 

Regardless of the language in which the computer program 
is written, it can be expected to contain some errors. Care- 
ful planning and attention to detail will go a long way to=- 
ward minimizing these errors, but a computer code is an ex~ 
tremely complicated structure, and the specifications for 
computer programs are very exacting. The omission of one ine 
struction or even one decimal point may result in the entire 
program being rejected by the computer. 

It is probably easier to test the program if it has been 
constructed in independent sections, and it was partly toward 
this end that the model was divided into functional components. 
Fach of these sections provides a particularly convenient unit 


when testing the program because of the relation between input 
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data and individual functional components. If the coding is 
tested in blocks, input data which would normally be provided 
by other sections of the program is usually introduced by a 
check-out control program. This may consist of a few cards 
containing the necessary data or might incorporate previously 
tested subroutines. In this way, data which is input to the 
subsection being tested can be precisely controlled. lLocal- 
izing an error in @ program section which uses self generated 
input data is complicated by the fact that both the program 
logic and input data are possible sources of the error. 

Although coding the program in subsections is probably 
the most significant single aid to program testing, it is not 
a complete panacea. Though each section performs properly 
over its designed range of inputs, this does not insure that 
all of these components will interact with one another as the 
designer intended. One method of testing the completed pro= 
gram is to compare a limited number of runs with hand calcu- 
lated results. This will generally involve checking the 
program at certain critical points. The procedure used in 
testing AHS-l1 was to have the values of key variables printed 
at certain points throughout the program. In this way the 
output from one operation and the input to the next link in 
the chain were checked simultaneously. These were then com= 
pared with hand calculations. 

The above procedure is a tedious and time consuming task 
for even a relatively uncomplicated program, but it may be 
the only way the simulation designer can be reasonably con- 


fident of the correctness of the program. Even when such a 
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testing procedure has been completed there is generally no 

way of knowing with one-hundred percent confidence that the 
program is free of errors. Such assurance would ordinarily 
involve testing every possible combination of input values, 
a virtually impossible task. 

5-4 Auxiliary program components. 

The questions of what to include as input and output data 
confronted the simulation designer throughout the model de- 
slgn phase. The question of how to get this information in 
and out of the computer must also be answerede With the cap= 
acity to process large amounts of data inherent in a digital 
computer, data processing is not a trivial consideration. 

The input data format will often be influenced to a cone 
Siderable extent by the size of the program. Space limita= 
tions may necessitate reading data in only as it can be used. 
If space is not a problem, time may be saved by reading as 
much data at one time as is practicable and storing it intern- 
ally. 

Space and time considerations will affect the particular 
form of input data as well. In a small simulation, data may 
be introduced in a form which is more convenient to the user 
and any necessary conversions accomplished internally. Large 
simulations might require considerable pre-processing of data 
to conserve room for calculations which can be conducted only 
in conjunction with the operating section of the progran. 

Errors in the input data are also an important consideraq= 
tion. Data input to AHS=-1l is printed out as a check on what 


values were actually used by the simulation, but an important 
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section of the program was omitted because of time limitations. 
This omission is a routine for checking the input quantities 
to see that they fall within acceptable limits. This would 
not eliminate mistakes entirely but could reduce lost comput- 
ing time due to misplaced decimal points and other gross errors 
in the input cata. 
The output section of the program will be influenced to 
& considerable extent by the simulation objectives. In addi- 
tion to selecting outputs which will yield the maximum informa- 
tion, the designer should specify an output format which is 
compatible with the particular analysis techniques to be 
employed. For example when space permits, computation of 
sample means and variances might save the analyst a great deal 
of time. Similarly, percentages or ratios of numbers are often 
more meaningful than the numbers themselves. 
5.5 Testing the simulation. 
In designing the model, factors were included because 
they were known, or in some cases only suspected, to affect 
the system in @ way which was important to the study purpose. 
Before the simulation is ready for use as an analytic tool, 
the designer should determine both quantitatively and qualita- 
tively how these factors affect simulation results. This will 
require that some testing program be conducted and representae= 
tive levels of all or part of the factors examined. 
Requirements of a testing program. The testing program 
may include a comparison of simulation results with observa- 
tions of the system itself, or with experimental results ob-= 
talned by other analytic methods. It should include some an- 


alysis of the sensitivity of measures of effectiveness to 
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changes in factor levels. Whatever the analytic ‘echniques 
employed, the following appear to be fairly general objectives 
of the tests. 


(1) Verify that simulation response to well known factors 
is in general agreement with observation. 


(2) Determine the sensitivity of response to all factors. 


(3) Determine the number of replications necessary to 
attain certain levels of statistical validity. 


(4) Determine which factors interact and the extent of 
such interaction. 


Conducting the testing program. The number of runs neces- 
sary to attain the desired accuracy should be one of the first 
considerations during simulation testing. In testing the simu- 
lation, as well as in conducting subsequent experiments, the 
analyst is usually guided by two major considerations; how to 
gain as much useful information as possible from the simulation, 
and still minimize the expenditure of time and resources. Test= 
ing should also provide some general guidelines as to the number 
of replications necessary to achieve certain levels of confidence 
in the simulation results. These may then be applied during both 
employment and testing of the simulation. Testing the simula- 
tion using too few runs may essentially invalidate the results 
of the tests. 

When comparing simulation results with observations of the 
system, it 1s well to consider the source of such observations. 
In analyzing military operations the results of peace time ex-= 
ercises may be the nearest the analyst can attain to an actual 
observation. Field maneuvers and fleet exercises may at times 


approach reality, but the fact remains that they are another 
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form of simulation. As with real battles, an exercise may be 
conducted only a few times. Consequently the analyst must con- 
sider the dangers associated with drawing inferences from small 
samples. The generality of the model must also be considered 
when ccmparing simulation response with observations from the 
real world. agreement between simulation results and one parti- 
cular observation may reflect that the designer has only simu- 
lated that particular observation. 

During the initial testing it may be noted that certain 
combinations of factor levels tend to saturate the simulation 
Output. That is, some factors because of their greater effect 
on the simulation response, may tend to completely determine 
the results when assigned values near the extremes of their 
ranges. Knowledge of these saturation levels should be use- 
ful when corducting sensitivity studies on the remaining input 
parameters and variables. 

Since the simulation was designed as a factorial experi- 
ment from the start it is expected that this method of analysis 
will play an important part in the sensitivity testing. Factor- 
lal experiments are particularly adaptable to sensitivity anal- 
ysis as well as to the determination of factor interactions. [10,26| 

Re Je Matteis and W. C. Sualer have described a testing 
program employed in sensitivity testing of the Carmonette Model. 
[26| The testing program was conducted in two phases; The first 
phase consisted of examining the more important factors and gen- 
erating some estimate of effects and variances. In this phase, 


factors which exhibited saturation effects were treated 
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individually. The remaining factors were then examined in a 
complete factorial design. In the second phase, knowledge 
gained during the initial testing was used in designing a frac- 
tional factorial experiment to test the overall model for inter- 
actions. 

5.6 Documentation. 

Documentation can very easily be the defining line between 
a partially completed computer simulation and a useful analytic 
tool. The importance of assumptions has been emphasized through- 
out the preceding discussion of model design. Unless these ace- 
sumptions are written down, however, they may not be considered 
during the analysis. This is further complicated by changes in 
personnel or by physical separation of designer and ultimate 
user. 

Properly documenting program logic and computer code is no 
less important than completely describing the model. The real 
world is dynamic and the simulation must often reflect change. 
Attempting to modify a program without complete flow charts and 
a dictionary of variable names, as well as the major model as- 
sumptions, can be a formidable task. In some cases it may be 
simpler to rewrite the coding than to modify an existing program. 

Standards of format and adequate content vary with each or- 
ganization engaged in simulation design, and standardization 
among simulation groups is virtually non-existant. The descrip- 
tion of AHS-1 reflects the author's own idea of adequate docu- 
mentation, with one exception; the simulation has not been 
thoroughly tested and no mention is made of tests results. In 
summary, the following list is considered to contain the minimum 
requirements for documentation of a computer simulation. 
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(1) 


(ay) 


(3) 


(4) 


©) 


(6) 


A statement of all major assumptions pertaining to 
the model with appropriate references, and a descrip- 
tion of the model. 


A description of the logic. This should include gen- 
eral flow charts. 


Complete rules for input data preparation, and the 
allowable range of input values. 


Statement of measures of effectiveness and a descrip- 
tion of the output data. 


A list of all variable names used in the program with 
their definitions. 


A description of the statistical test procedures used 
and the results of these tests. 
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CHAPTER VI 
AHS-1 4&4 COMPUTER ‘SIMULETION 


AHS-1 is an event store computer simulation of an antieg 
submarine warfare helicopter small area search. In accordance 
with the definitions of Chapter I, it may properly be described 
as an analytic computer war game. A maximum of ten dipping 
sonar equipped helicopters and one evading submarine are the 
principles in the game. 

The simulation is intended as a tool for the statistical 
analysis of the comparative effectiveness of different helicopter 
search tactics in similar tactical and physical environments. 
The model for AHS-1l and the underlying assumptions concerning 
sonar parameters and submarine capabilities are described in 
detail in the succeeding sections. 

6.1 fJactical situation and play of the game. 

A play of the game begins with the submarine at a known 
datum and a flight of sonar equipped helicopters enroute to or 
in the vicinity of the datum. The submarine leaves the surface 
at game time zero and departs the datum on a course and speed 
Which is unknown to the helicopters. The helicopters arrive 
at datum at some time determined by the game user and proceed 
to dip stations in accordance with the assigned search plan. 

The submarine leaves datum on a course, speed and depth 
selected in one of four ways according to the game inputs. 
Fach of these operating modes reflects a different degree of 


submarine randomness, ranging from a completely predetermined 
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track to random selection of course, depth and speed within 
restrictions imposed by the game planner. 

The helicopters proceed to their designated dip stations 
at datum time plus time late and commence the search. At each 
dip station as many as five sonar sweeps at various depths can 
be executed. Helicopters are assumed to be equipped with scan- 
ning type variable depth active sonar. The time to complete 
each sonar dip is determined by input time delays representing 
the following: (1) The time initially required to establish 
a hover and to lower the transducer, (2) the times necessary 
to change the transducer depth between sweeps, (3) the elapsed 
time from completion of the last sonar sweep until the helicopter 
has transitioned from hovering to forward flight, and (4) time 
spent in the prosecution of non-submarine contacts if applic=~ 
able. The dip cycle time is the sum of dip time and time ene= 
route between dips. Except for delays resulting from non=<sub- 
marine contacts, dip time is the same for each helicopter. All 
helicopters fly at the same airspeed between dips. 

At the beginning of each sonar sweep the submarine position 
1s determined and a target range computed. The detection area 
is doughnut shaped, centered about the helicopter'’s position. 
Detection will occur if the submarine is within the annulus 
determined by the maximum detection circle and a smaller circle 
Within which detection cannot occur because the helicopter is too 
Close to the target. This minimum detection range is not de= 
pendent on sonar conditions and may be input as zero if desired. 


The maximum detection range is computed for each sweep and each 
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helicopter at the beginning of the play and recomputed whenever | 
the submarine changes depth. Detection ranges are based on a 
randomly selected figure of merit for each unit and the prevaile 
ing sonar conditions. Once selected, individual helicopter 
figure of merits are constant throughout each play. 

The scheme for computing detection ranges was designed by 
R. Lo. Klinkner of the Applied Physics Laboratory [ 17] and is 
considered by the author of this thesis to be an important im= 
provement over the deterministic “cookie-cutter” detection range 
employed by many sonar simulations. The actual detection range 
for each helicopter depends on sonar frequency, sea state, layer 
depth, temperature in the layer, target depth, transducer depth, 
and the quality of the sonar operator-equipment combination. 

The latter is assumed to vary among helicopters in a prescribed 
random manner. 

The play proceeds until all helicopters have completed 
their last dip or a detection occurs. At the completion of the 
play the running tally of detections is brought up to date and 
anew play commenced. This sequence continues until the desired 
number of plays have been completed. Input parameters and vari- 
ables may then be changed and another sequence repeated. A 
series of plays with any given set of input data will be refer= 
red to as @ game run. 


6.2 Submarine movement. 





The target submarine maneuvers in 4& three dimensional playe- 
ing area. Horizontai movement is relative to an X,Y coordinate 


System, with yards the basic unit of distance. Depths are 
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measured in feet to the nearest foot. Courses are to the near=- 
est degree and speeds to the nearest knot. Position and speed 
vectors are computed using a standard polar coordinate system; 
however, courses and speeds are adjusted so that the Y axis cor- 
responds to north and all courses are measured clockwise from it. 

Four different modes of submarine operation are available to 
the user of AHS-1 as follows; 

Option (1). Submarine track is predetermined by the user. 
This mode might be used when studying the effect of evasive sub- 
marine maneuvers on a particular helicopter search plan. If the 
probability of delay due to non-submarine contacts (see Section 
6.4) is input as zero, the location of the helicopters will be 
known at all times to the game planner. Any level of intelli- 
gence concerning helicopter movement may then be attributed to 
the submarine. Similarly, dip times may be made essentially 
random by assigning positive probability to non-submarine con= 
tacts. In this way the submarine may be placed in the situation 
of knowing where the helicopters are at any one time but not 
knowing when they will move, or where they will jump to. 

Option (2). Speed, depth and the bearings of each track leg 
relative to the first course are predetermined as in Option (1). 
However, at the beginning of each play a pseudo-random number 
1s selected from the interval [o, 360) degrees and added to each 
course. Election of this mode is equivalent to assuming a come 
plete lack of knowledge, by the submarine, concerning helicopter 
movement. AS opposed to a completely random submarine (Option3), 


however, the capability of changing course, depth and/or speed 
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during the play is retained. 

Randomizing the submarine track serves to preclude accident- 
al biasing of game results by the user. This bias may be in-= 
troduced by a particularly fortuitous selection of the sub- 
marine track with respect to the helicopter positions. A site 
uation particularly suited to this mode is that of a submarine 
departing datum on a straight course but decreasing speed at the 
end of preset time intervals. This tactic might be used by a 
submarine commander who desires to open the range to datum as 
expenditiously as possible without exhausting the submarine 
batteries. 

Option (3). Submarine course, speed and depth are uniformly 
distributed random variables. Course is distributed between zero 
and 360 degrees. Depth and speed range between upper and lower 
limits chosen by the game user. Course, speed, and depth are de= 
termined by generating three pseudo-random numbers at the begin= 
ning of each play, and remain the same until the termination of 
that play. 

Ortion (4), Course and speed are chosen randomly as in 
Option (3) but depth is uniquely determined by the chosen speed 
and the game inputs. After selecting a speed, the minimum depth 
is taken as the shallowest depth at which the submarine can oper- 
ate without cavitating. It should be noted here that cavitation 
does not directly affect the possibility of detection in the 
game, as the helicopters do not conduct passive sonar operations. 
This submarine mode is introduced as a means of restrictins the 


target to depths which are realistic for the speed used. ‘arget 
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depth is a parameter which is considered in determining active 
sonar range. 
6.3 Helicopter movement. 

Helicopter movement consists of jumps between sonar dip 
stations. All helicopters jump at the same speed, and as the 
game is now programmed dip stations are deterministic. Provisions 
have been made in the model and accompanying computer code for 
the introduction of random bearing and course errors. However, 
determination of the distribution of such errors would require 
the services of fleet units may’ available to the writer. With 
the increasing cond b thedh tee Wt helicopter navigation equipment 
the omission of these errors is not considered to be critically 
detrimental to the purpose of this simulation. The provision 
for including navigation errors is intended primarily as an area 
for further study by the interested reader. 

Helicopter movement over the playing area is unrestricted 
in azimuth, and dip stations are determined by the game user. 
The game is primarily intended to simulate close area search 
plans, and helicopter dip stations are computed relative to 
datum. When simulating screening or other support operations 
the datum coordinates may represent an aircraft carrier or other 
helicopter take-off point. The first submarine maneuver may 
then be used to move the submarine from datum to the desired 
starting point. This maneuver can be executed at a high enough 
speed so that it will be essentially instantaneous. 

6.4 Submarine detection. 
rach helicopter has the capability of making up to five sonar 
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sweeps at each dip station. Transducer depth and the time re- 
quired to adjust transducer depth and conduct a search are in- 
put by the user and may be different for each sweep. The number 
of sweeps, sensor depthandtimetocomplete corresponding sonar 
sweeps are the same for all helicopters. 

No provision is made for submarine motion during a sonar 
sweep. This is of very little consequence if scanning type sonar 
1s simulated, but could have some effect on the outcome if search- 
light sonar is being used. In the case of searchlight type sonar, 
the time to train the transducer may be included in each sonar 
sweep time; however, range to the submarine will only be deter- 
mined at the beginning of each sweep, regardless of the time re- 
quired to step train the transducer through a complete sweep. 

At the beginning of a sonar sweep the range to the submarine 


1s computed according to 


i” ff zl a  ee 
an 2 _ 2 
Rs =)/ (Xp, Ral” + (po Kee) (Gu) 
Where xX and Yh are the dip coordinates of the helicopter, 


and x, and We are the submarine position coordinates. Detec=- 
tion occurs if SR = Rs = Rd where Rd is the detection range 
for the particular helicopter at the appropriate transducer 
depth, Rs is the target range, and SR is the minimum range at 
Which detection can occur. SR may be input as zero if desired. 

The course, speed and depth of the submarine at the begin- 
ning of a dip are used in computing target range. If the target 
is scheduled to maneuver during any sweep this information is 


noted and at the completion of the sweep in question the dip is 
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discontinued. The submarine maneuver is then executed and the 
dip completed using the correct submarine track information. 

In addition to sonar sweep times and delays representing 
the time required to lower and retrieve the sonar transducer, 
one other factor contributes to sonar dip time. This is a random 
variable representing time engaged in pursuing non-submarine con- 
tacts. 

For the purposes of this game the time required to classify 
actual submarine targets is not considered pertinent. If a sub- 
marine has been detected, the search plan has been effective and 
classifying the target is another problem. It is conceivable, 
however, that the effectiveness of different search plans may be 
more or less sensitive to variation in dip times among helicopters. 
One of the primary factors contributing to differences in dip 
time is the classification of false contacts. In the model it 
has been assumed that non-submarine targets such as fish or 
whales are uniformly distributed throughout the ocean. This im= 
plies that one sonar operator has about the same chance of con-= 
tacting such a target as any other sonarman. Once a non-submarine 
contact has been generated, the time required to classify it 
as non-submarine will vary among sonar operators. These two 
assumptions have been incorporated into the model in the follow- 
ing way. 

The actual delay for each helicopter is dependent on two 
input parameters, the probability of obtaining a false contact 


(P,,) and the maximum time any non-submarine contact will be 
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prosecuted (Tmax). A random number from the interval lo, 1] is 
first generated to determine whether the helicopter gained a 
false contact during the dip in question. If this random nun- 
ber is less than or equal to Pro a non-submarine contact was 
acquired and a delay time must be added to dip time. The length 
of time to be added is determined by 


T = — . Tmax (i6m2)) 


1 iG 


Where T is the actual delay and RN is the previously generated 
random number. 
6.5 Sonar detection ranges. 

Detection ranges for each helicopter are computed at the 
beginning of every play and whenever the submarine changes depth 
during the play. As detection range depends on transducer depth, 
every helicopter will have associated with it as many different 
detection ranges as the number of sonar sweeps per dip. The 
actual computation of these detection ranges is a problem in 
acoustics rather than simulation and will not be considered in 
detail. The interested reader is referred to Fundamentals of 
SONAR, by J. W. Horton [19] for further information on this sub- 
ject. 


The basic equation for active sonar is; 
Ep =" Fi Ss - 2EL (Gna) 


where EE is the echo excess (signal level relative to that re= 


quired for a 50% probability of detection) 


PM = Sonar figure of merit 
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TS = Target strength 

PL = One-way propagation loss 
All quantities in equation (6.3) are expressed in decibels. 
Figure of merit and target strength are input variables. Figure 
of merit is a measure of the quality of the sonar operator- 
equipment combination independent of water conditions. Target 
strength is a characteristic of the submarine size and shape, the 
external surface of the submarine hull, and target aspect. An 
average value of target strength is used in the model. An echo 
excess of zero corresponds to a 50% probability of detection 
and this value is used in the computation. 

In computing sonar ranges it is assumed that sonar figure 
of merit varies among helicopters according to some known prob= 
abllity distribution. It is also assumed that figure of merit 
does not vary appreciably between dips for the same helicopter. 
Therefore the figure of merit is computed for each helicopter 
Only once for each play of the game. In order that this thesis 
be unclassified no attempt has been made to duplicate the actual 
distribution of helicopter figure of merit. For purposes of 
AHS-1 a Normal (gaussian) distribution is assumed with mean and 
variance as input parameters. 


Rewriting equation (6.3) with echo excess as zero: 


Propation loss (PL) combined with transducer depth and the re- 
maining sonar parameters, is used by the routine which computes 


detection ranges for each individual helicopter. MThe other 
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input variables entering into detection range are target depth, 
sonar frequency, layer depth, temperature in the layer and sea 
state. 
6.6 Description of program subroutines. 

The following brief descriptions of the individual sub- 
routines are intended to show how the model described above 
has been implemented. A secondary purpose of these descriptions 
is to indicate auxiliary functions performed by each event sub- 
routine. Rules for input data card preparation are included 
under Subroutine PRINT. Flow charts showing the logical struce- 
ture of each subroutine are included at the end of the individual 
descriptions. the Following flow chart symbols are used through- 
Out. 

Flow chart symbols. The meaning of each symbol is determined 
by shape and letters within the symbol. Differences in size re= 


flect only space consideration. 


Beginning point in a subroutine. 


- Computation to be carried out within 
the subroutine, 


- Arrows show direction of flow, 


~ Computation to be carried out by a sube- 
routine other than the one in which the 
symbol appears. In the detailed flow 
charts the information needed by the 
routine being called is written inside 
the block. 
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RETURN 


STOP 


-) 


} 
O 
MN 


(| i ff} Bo © O® 


Decision point in the logic. The route 
to be followed is indicated by the YES 
or NO. 


A logical jump is to be made to some 
point within the subroutine. 


Beginning point of a subsection of logic 
Within a subroutine. A logical jump will 
always go to such a point. The numeral 
indicates page number. 


Continuation of the logical flow to the 
next page. This symbol would appear at 
the bottom of page one and at the top of 
page two. 


and of event subroutine. Program is to 
process the next event. 


End of an event subroutine when for any 
reason the normal program logic is to be 
interrupted. 


Logical end of any subroutine which is 
not an event routine. Indicates a jump 
back to the section of the program that 
called the subroutine into action. 


Instruction to stop the game. Indicates 
that all the data has been processed or 
a mistake was made in the input informa- 
tion. 


Data is to be read into the machine. In 
the detailed flow charts the top numeral 
indicates the tape unit to read from and 
the characters in parentheses indicate the 
format to be used. NAME is the variable 
name to be read. 
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- Write output on tape unit six. The nunm- 
eral at the bottom is a code indicating 
the format. 


SNE -~ Store New Event. Type of event and 
event time will appear within the 
block. 


Symbols used within blocks in detailed charts. 





RN - Random Number 
- Multiplication 
/ - Division 
4H - Exponentiation 
A+B = C - Add A to B and store the answer in the 


variable name C. 

List of variable names. Following is a list of the variable 
names With their definitions. The names are arranged according 
to the subroutine in which they were originally defined. This 
also corresponds to the order of the COMMON statements. Input 
data names are designated by an asterisk. Names representing 
data to be printed are identified by $ or # symbols. Dimensioned 
variables are indicated by parentheses following the name. 


(1) Utility Words 


I™ Used for temporary storage or con-= 
Ti. venience throughout the program. 

Te 

ONG Helicopter dip counters. I2 is set 
Ie equal to the total number of dips 


(number of helicopters times num-= 
ber of dips per helicopter) at the 
beginning of each play. IDIP is 
incremented each time a helicopter 
dips. If a detection does not occur 
the play is over when IDIP = [I2. 
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Je 


Counter for number of samples. Used 
in computing average number of de- 
tections per sample. 


(2) Main Control Program 


#NRS AM 


*NSSIZE (I) 


$NHSDET (I) 


$TMIN 


¢TMAX 


¢BART 


$VART 


#MIND 


HM AXD 


#BARD 


Number of samples in each simula- 
tien rum. 


Number of plays constituting the 


‘len sample. 


Counter for keeping track of the 
number of detections made during 
one sample by the pth helicopter. 


Minimum time in minutes from the time 
the helicopters arrive at datum until 
a detection occurred for current 
sample. 


Maximum time required for a detece= 
tion during any sample. 


Average time elapsed between time 
helicopters arrive at datum and 
detection when detection occurs. 


Sample variance of detection time 
when detection occurs. 


Minimum number of plays resulting in 
detections during a series of runs. 


Maximum number of plays resulting in 
detections during a series of runs. 


Average number of plays resulting 
in detections per sample during a 
series of runs. 


Sprinted after each sample is completed. 


ft 


Printed at the end of runs designated by game user. 
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# VARD 


$ PERC 


* DAT 


* XDAT 
* YDAT 
$ NDET 


NDTEMP 


PL(I) 


* IRNC 


Sample variance of number of detections per sample. 


Percentage of plays resulting in a detection dur- 
ing a sample. 


Time at which helicopters arrive at datum and 
begin search. 


X and Y coordinates of datum. 


Counter which is set to zero at beginning of each 
sample and incremented each time a detection occurs. 


Temporary counter which is set equal to NDET at the 
beginning of each play. Comparison of NDET and 
NDTEMP signals whether the play ended because a 
detection occurred or because the helicopters 
completed all dips. 


One way sonar propagation loss of the ye helicopter. 
Computed at the beginning of each play. 


Number of pseudo-random numbers to be generated 
and disposed of at the beginning of a simulation run. 


(3) Subroutines SNE and TNE 


The following three words are used throughout the program 
to transfer information to or from the Event Store Table. 


TIMET 
NREVT 
NRUNT 
NTNE 

ed, 


NREV 


1 
I 
NRUN(I) 


Time at which an event is to be stored or executed. 
Number of the event to be stored or executed. 
Number of the unit affected by the event. 


Counter which keeps track of the location of the 
last event in the table. 


These three arrays comprise the event store table, 
and represent time of execution, number of the 
event subroutine to be called, and the number of 
the unit involved in the event. 
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(4) 


s¢ 


+ 


KT 
NRANSS 


NRANC 


oTIME 


SMNEXT 


It 
VAS 
VY5 
AS 
YS 
bi 


Subroutine SM= (Submarine Maneuver Event) 


Submarine maneuver counter. 


Signals whether submarine maneuvers are predeter- 
mined or random. 


Signals whetner the first course is randomly se- 
lected when submarine maneuvers are predetermined. 


Time the submarine jtast maneuvered. 


Time at which the submarine is scheduled to make 
the next maneuver. 


Depth at STIME. 
a at time STIME. 


Vy at time STIR. 


X and Y coordinates at time STIME. 


Temporary storage ror random numbers which are to 
be added to submarine courses. 


Submarine maneuver table. 


$¢ 


3t 


3t 


4t 


SMTIME(I) 
DEPTH(I) 
S6USs(1) 
SSPD(I) 
NRM AN 


h 
t submarine maneuver. 


Time to execute the I 
Depth (feet) at time SMTIME(I). 

Course (degrees) at time SMTIME(I). 
Submarine speed (knots) at time SMTIME(I). 


Number of times the submarine is to maneuver 
when maneuvers are predetermined. 


Input parameters when submarine maneuvers are randon. 


+ 
+% 


a 


* 


NHISPD 
LOSPD 


MAXSD 
MINSD 


Upper and lower Limits (knots) of submarine speed. 


Upper and lower limits (feetjof submarine depth. 
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Cavitation speed versus depth table. 


%* NCAVS(I) 


# NC AVD(I) 


*# NC AV 


Cavitation speed in knots for the 1th point 
on the curve. 


Minimum depth corresponding to NCAVS(I). 


Number of entries in the cavitation table. 


(5) Subroutine HSM (Helicopter Maneuver Event). 


IDIPN(I) 


NRDIPS 


HTE(I,d) 


HBRG(I,J) 


XH(1) 
WH (1) 


NRHS 


HSPD 


HSPDT 


TDD 


ERB 


ERT 
XHT(TI) 
THe( £) 


NTE 


Counter which keeps track of the number of dips 
made by the Ith helicopter. 


Number of dips to be made by each helicopter. 
Jump time of [th helicopter enroute to the gth 


dip station. Input as yards and converted to 
minutes internally. 


Bearing relative to the preceding leg flown by 
the I&" helicopter in transiting to the J°&" dip 
station. Input as degrees and converted to 
radians internally. 


X and Y coordinates of the Ith helicopters last 
dip station. 


Number of helicopters. 


Speed in knots at which helicopters transit 
between dip stations. 


HSPD converted to yards per minute. 


Time (minutes) required to establish a hover 
and lower the transducer. 


Error in bearing added to HBRG(I,J) to de- 
termine actual flight path of helicopters. 


Distance error added to HTE(I,J). 

Table for temporary storage of helicopter X 
and Y coordinates. Used whenever ERB and ERT 
are not computed. 

Signals whether ERB and ERT are to be computed. 


Value is one if errors are to be computed, and 
zero otherwise. 


Ge 





(6) 


(7) 


DRNG(I,J) 
NRSW 
SWT(I) 


PD(T) 


TRD 
NOVER 
SHORTR 


NSWP(T) 
PRFO 


FC TM ax 


ICT 


NS 
> 
KF 
BP 
DL 
TS 


FOMMU 


FOMV AR 


Subroutine HDE (Detection Event) 


Detection range of the pth helicopter for the 
Jth sonar sweep. 


Number of scnar sweeps to be made by each 
helicopter at each dip station. 


Time (minutes) required to complete the ia 


sonar sweep. 
Transducer depth for the pth sonar sweep. 


Time (minutes) required to retrieve the trans- 
ducer at the end of each dip. 


Signals that a play has been completed. Value 
is one if play is over, zero otherwise. 


Range (yards) from helicopter within which 
target IS tod néer for d@tectifon to occur. 


Signals the first sweep number to be executed. 
Probability of detecting a non=-submarine contact. 


Maximum time for evaluation of non-submarine 
contacts. 


Actual time to evaluate non-submarine contact. 


Subroutine COMPRG 


Sea State. (beaufort scale). 

Temperature in the layer. (degrees Fahrenheit). 
Acoustic frequency of the helicopter sonar. 
Transducer depth. (feet). 

Layer depth. (feet). 

Target strength. (decibels). 


Average figure of merit of the helicopter sonar. 
(decibels). 


Variance of helicopter sonar figure of merit. 
(decibels). 


T) 





Main Program. The main AHS-1 program performs the functions 
of moderator and bookkeeper. This section of the program calls 
for more input data when it is needed and records the results of 
each play. The main program also sets up the initial submarine 
and helicopter maneuver events so that each replay of the game 
Will begin at the proper point in time and space. 

The following terms will be referred to throughout the de- 
scription of the progran. 

Play. A play of the game begins when the submarine leaves datum 

at game time zero and ends either when the submarine is detected 

by a helicopter or when all helicopters have completed their last 
dip. 

Sample. A sample consists of a series of plays. Sample size is 

determined by the user and is ordinarily based on what is consid- 
ered necessary for statistical validity of the results. 

Run- A run consists of one or more samples. 

At the completion of a run, any number of input parameters 
may be changed. Inputs do not vary between samples of the same 
run, additional samples representing replications of the same ex- 
periment. The only difference between two samples of the same 
run is the series of random numbers used. The random number gen= 
erator is reset to the first number in the series at the beginning 
of each run. 

The various functions of the main program should be apparent 
from the accompanying flow diagram once the reader is familiar 
With each individual subroutine. Only those parts which are not 


immediately clear will be mentioned here. 
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It is contemplated that the usual submarine operating mode 
will be a completely random submarine, option (3). The datum 
coordinates will normally be the origin of the X,Y, coordinate 
system. For this reason variable names associated with these 
inputs, referred to as optional variables in the flow chart, 
are preset prior to the first run. This obviates the necessity 
of reading these inputs in except when other than this standard 
operating mode is desired. 

The minimum number of detections and maximum time to detect 
are initially set at large numbers so they will always be replaced 
by the actual quantities once a detection occurs. If the sub- 
marine is not detected these qauantities will be printed as 1000 
detections and 600 minutes respectively and should be disregarded. 
If a sample size of 1000 or larger is used, the main program 
should be changed to reflect this by setting MIND equal toa 
number larger than the sample siZe. 

Bearing and course inputs to the program are measured clock- 
Wise, courses being relative to north,or 000. For computational 
convenience the zero radial is congruent with the X axis and all 
bearings are measured counter-clockwise. For purposes of the 
game itself this presents no difficulty since only relative dis- 
tances and bearings affect the simulation results. With no ad- 
justment the resulting game would be a mirror image of what one 
would plot from the inputs, and relative distances would be pre- 
served. In anticipation that the user might wish to have unit 
positions printed out at some time during the play, an adjust- 


ment has been made to all courses. By using the negative of 
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courses as input to the program and adding ninety degrees, the 
coordinate system has been transformed into a geographical map. 
Consequently submarine and helicopter positions will agree with 


the usual map representation. 
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MAIN PROGRAM (General ) 


| START 


SET OPTIONAL INPUT 


VARIABLES TO ZERO 





SET CONSTANTS FOR 
SONAR RANGE 
COMPUTATIONS 


SET FIRST ENTRY 
IN EVENT STORE 
TABLE TO ZERO 


INITIALIZE DETECTION 
STATISTICS 


SET GAME TIME TO ZERO (a2) 





Main Prog. 
a GF 6 







IS A NEW RUN REQUIRED 


READ IN DATA 
AND PRINT 
NEW RUN INPUTS 


INITIALIZE NEXT SUBMARINE MANEUVER 
TIME SIGNAL FOR RANDOM SELECTION 


OF SUBMARINE COURSE AND SPEED 





SET SUBMARINE COURSE ADJUSTMENT 
TO ADD 90 DEGREES TO ALL COURSES 


ARE RANDOM ERRORS 
TO BE ADDED TO YES 
HELICOPTER JUMP 
BEARINGS AND RANGES 


0 





COMPUTE AND TABULATE 


DIP STALIONS 
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TO 
B3 





Main Prog. (G) 
© 5 ef © 


(23) INCREMENT SAMPLE COUNTER 


RECORD AVERAGE FIGURE OF 
MERIT AND FOM VARIANCE 









RECORD TOTAL NUMBER 
OF DIPS AND SET TOTAL 
DIP COUNTER TO ZERO 


RESET AND CYCLE 


RANDOM NUMBER 
GENERATOR 


BEGIN NEW (3. 
SAMPLE 
INITIALIZE INDIVIDUAL 
HELO DIP COUNTERS 





(2 





C4 Main Prog. (G) 
4 on 6 





INITIALIZE MINIMUM 
MAXIMUM, MEAN AND 
VARIANCE FOR 
DETECTION TIME 






BEGIN | (oi 
NEW PLAY 


INITIALIZE SUBMARINE DEPYH 








SET UP FIRST 
SUBMARINE MANEUVER 
AND STORE IN EVENT 
STORE TABLE 











SBT UP Piet BrP 
FOR EACH HELICOPTER 

AND COMPUTE PROPAGATION 
LOSS FOR EACH HELO 


INITIALIZE END OF GAME FLAGS. 
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(4% Main Prog. (G) 
5S of 6 










TAK= FIRST EVENT, 
TNE WILL HAVE CONTROL 
UNTIL DETECTION OR 
COMPLETION OF LAST DIP 


NO WAS SUBMARINE DETECTED 


YES 


INCREASE INDIVIDUAL HELO 
DETECTION COUNTER 


COMPUTE SQUARE OF TIME TO 
DETECTION AND RECOMPUTE 
MINIMUM AND MAXIMUM 
TIME TO DETECTION 





IS THIS THE LAST NO 
PLAY OF THIS SAMPLE 


YES 
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COMPUTE MEAN AND VARIANCE 





OF TIME TO DETECTION 









COMPUTE PERCENTAGE 
OF PLAYS RESULTING 
IN SUBMARINE DETECTION 













COMPUTE SQUARE OF NUMBER 
OF DETECTIONS AND RECOMPUTE 
MINIMUM AND MAXIMUM OF 
DETECTIONS PER SAMPLE 


IS THIS THE LAST NO 
AMPLE FOR THIS RUN 


ED 


GD 
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Main Prog. (G) 
6 of 6 


Q> 





Subroutine SNE. Subroutine SNE (Store New Event) enters 
future events in chronological order in the event store table. 
The event store table is a list, ordered according to time, of 
actions which are to take place in the course of a play of the 
game. 

When an event is to be stored, SNE is entered with three 
items of information; the time at which the event is to occur, 
the type of event to be executed, and the number of the unit in- 
volved. This information constitutes one event word. 

Should two or more events be scheduled to occur at the same 
time there is no designated ordering according to event type. 
Two or more such events will take place in the order in which 


they were processed by SNE. 
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Subroutine SNE (General ) 








COMPARE EXECUTION TIME 
OF EVENT TO BE STORED 
WITH EVENT TIMES IN 
EVENT STORE TABLE 









INSERT EVENT INFORMATION 
(TIME, UNIT NUMBER, EVENT 
NUMBER) IN PROPER 
CHRONOLOGICAL LOCATION 
IN EVENT STORE TABLE 



















INCREASE EVENT 
COUNTER TO REFLECT 
ADDITION OF EVENT 
TO EVENT STORE TABLE 





Subroutine SNE (Detailed) 


I = NINE 
NTNE = N 


NO 
TIMET < TIME(TI) 














TIME(I) 
NRUN (I) 
NREV(I) 


TIME(I+1) 
NRUN(I+1) 
NREV(I+1) 









TIME(I+1) 
~ a Sed 
T+1) 


NREV 
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Subroutine TNE. Subroutine TNE (Take Next Event) examines 
the event type of the earliest entry in the event store table 
and calls the appropriate event subroutine into action. Once a 
play has commenced, the program is controlled by TNE until the 
end of the play. The end of a play is signaled by the detection 
event subroutine. When this signal is received, TNE transfers 
control back to the main program. 

When an event is executed the event store entry referring to 
it is removed from the table. An alternative method would be to 
leave the entries in the table and keep track of the earliest un- 
executed event. The method used is somewhat more time consuming 


but results in a considerable saving of storage space. 
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Subroutine TINE (Zeneral) 











HAS A 
DETEC TION 
OCCURRED 


RETURN 
CONTROL 


TO MAIN 
PROGRAM 





NO 


STORE EARLIEST 
EVENT WORD IN 
TEMPORARY STORAGE 


SHIFT STORED EVENT 
WORDS DOWN ONE CELL 
IN EVENT STORE TABLE 


CALL NEXT EVENT 
SUBROUTINE CORRESPONDING 
TO EVENT WORD IN 
TEMPORARY STORAGE 





87 





Subroutine TNE (Detailed) 


LES 


NOVER 2 0 (a) 


NO 














TIME(2) => TIMET 
NRUN(2) <>» NRUNT 
NREV(2) <>» NREVT 










TIME(I) <> TIME(I-1 
NRUN(I) <> NRUN(I-1) 
NREV(I) <> NREV(I-1) 





Subroutine SME. The first time this subroutine (Sub= 
marine Maneuver Event) is called on during a play, the subd- 
marine operating mode is determined from the information in- 
put by the game planner. Once this determination has been 
made, course, speed and depth are determined randomly or read 
directly from the input information. Game time and the sub- 
marine position are noted and xX and Y components of velo- 
city are computed. From these four pieces of information the 
submarine's coordinate position can be determined at any later 
time. 

Additional functions performed by Subroutine SME are 
as follows: 

(1) On the first submarine maneuver or any subsequent 
maneuver involving a depth change, helicopter de- 
tection ranges are computed by calling Subroutine 
COMPRG. 

(2) If another submarine maneuver is to take place, the 
time the maneuver is to occur is noted and subroutine 
SNE is instructed to place a submarine maneuver event 
in the event store table. 

(3) The time of the next submarine maneuver is recorded 
for use by the detection event subroutine. This in- 
formation is used in determining whether the sub- 
marine maneuvered during any helicopter dip. In the 
event of a submarine maneuver while a helicopter is 


dipping, submarine track and depth information must 
be updated before the dip is completed. 
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Subroutine SME — Submarine Maneuver Event (General ) 


ENTER 


NO IS THIS THE 
FIRST MANEUVER 


YES 
ens ARE SUBMARINE \<E8_-Go 2 
IS THIS A DEPTH COURSES RANDOM B2 
CHANGE ONLY 
IS PIRST COURSE 


RANDOM 











COMPUTE RN 
BETWEEN 
0-360 


ADD RN TO 
ALL COURSES 












COMPUTE Nw X, 
TV Vy and 


RECORD TIME 





COMPUTE 
io THIS & YES HELICOPTER 


DEPTH CHANGE DETES TION 
RANGES 









ONE 
EVENT 1 AT 

NEXT MANEUVER 
TIME 


NO 





NO 











IS THis THE 
LAST MANEUVER 
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SOMPUTR RN 


BETWEEN 
0-360 





SUBMARINE 
OURSE = Ri 


IS SUBMARINE 
NON-C AVITATING 


NO 


COMPUTE RN 
BE TW TEN 
SPEED LIMITS 


VW. 
SPEED = MINIMUM 
SFrEeD * RN 





COMPUTE RH 
BETWEEN 
DEPTH LIMITS 











V 
DEPTH = MINIMUM 
DEPTH + RN 






YES <0, 


OF 


COMPUTE VELOCITY 
(D2) COMPONENTS 
| Ve and om, 


SME (G) 
2 08 @ 





COMPUTE RN 
Be TWEEN 
SPEED LIMITS 






MINIMUM 


DD + RN 


DETERMINE 
MINIMUM 
DEPTH FOR 
NO CAVITATION 





COMPUTE RN 
BETWEEN 
DEPTH LIMIT 









DEPTH = MINIMUM 
DEPTH + RN 










JOMPUTE 
HELICOPTER 


DET TION 
RANGES 








Subroutine SME -Submarine Maneuver tvent (Detailed) 


ENTER 


- (01) 360 —> aNd 

<> CALL RAND 
YES 

rity 2 ( NRANSS=0 RN/57.29578 > Tl | 


NO 
NO 
YES 
U iO NRANC=O 
== NHISPD-LOSPD 
{+l —~> RNO 
CALL RAND 
60 -—> RNC 
CALL RAND (LOSPD+RN ) #33 .766 
— T2 








MAXSD-MINSD 






—_ TT 


+1 —_ RNC 
CALL RAND 


MINSD+RN —> DT 





oe 








NO 


‘cos 
SSPD(KT)#33.7666666 ->SPD | 


(=SOUS(KT)+TT)/57.29578 -> cUS 
sae aerate cota ae —> xs 
YS+(SMTIME(KT)=STIME)xVYS —> YS 
SPDx cos(CUS)—> VXS 
SPD# sin(CUS)-—> VYS 
SMTIME(KT)—> STIME 










DEPTH(KT) = Dt 


DEPTH(KT)<—> DT 






CALL 
-COMPRG 


\ 







KT+l —> KT 
SMTIME(KT)—> TIMET 
1—> NREVT 

CALL SNE 


a2 





SME (Dp) 
> ‘ote e 







NC AVS(NC AV )-NC AVS(1) 
+1 => #£=xRNC 
CALL RAND 









Toxncos(Tl)—> VXS 
To#sin(Tl)—> vVyYs 





| NOAVS(1)+RN —> NSS 
NSS#33.7666 =—> T2 


C ALL 
COMPRG 
l1-_- I 
SS — 
NSS < NCAVS(1) 


ihn 


NO AVD(NC AV)=NC AVD(TI) 
+1 -—> RNC 
CALL RAND 


NCAVD(I)+RN —> DT 
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Subroutine HSM. Subroutine HSM (Helicopter Maneuver Event) 
updates each helicopter's dip station coordinates and records 
this information for use by the detection event subroutine. As 
the program is presently written, dip coordinates are computed 
by the main program and stored in a table. This table is avail- 
able to Subroutine HSM whenever a helicopter is to be maneuvered. 
This requires that the computations be made only once for each 
game run. Should the necessary code for computing distance and 
bearing errors be written, new dip coordinates will be computed 
by Subroutine HSM every play. 

Additional functions of this subroutine areas follows: 

(1) Records the number of the last dip station occupied by 
each helicopter. After the first dip, the order of 
maneuvering is determined by dip cycle time rather 
than helicopter unit number. Therefore an individual 
dip counter must be maintained for each helicopter. 

(2) Stores the next detection event for the helicopter 
which has maneuvered. The succeeding detection event 
Will occur at the time the maneuver is completed plus 
the time which elapses while the transducer is being 
lowered into the water. The latter time includes normal 


delay time required to transition from forward to hover- 
ing flight and is an input to the simulation. 
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Subroutine HSM - Helicopter Maneuver Event 


DETERMINE NUMBER 
OF HELICOPTER 
TO BE MOVED 


INCREMENT 





DIP COUNTER 


ARE ERRORS IN 
COURSE AND DISTANCE 
TO BE COMPUTED 


YES 


COMPUTE COURSE 


AND DISTANCE 
ERRORS 





JOMPUTE NEW X AND 
Y COORDINATES 
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(General ) 





NO READ NEW X AND 
Y COORDINATES 


FROM TABLE 









HSM (G) 
o “ent 2 


COMPUTE 
TIME ENROUTE 





| ADD TIME ENROUTE 


AND TIME TO 
LOWER TRANSDUCER 
TO GAME TIME 


SNE 
DETECTION EVENT 
FOR THIS HELO TO 


\ OCCUR &T TIME 
\ COMPUTED ABOVE 


INE 


e) fi 





Subroutine HSM — Helicopter Maneuver Event (Detailed) 


ENTER 


NRUNT — I 
IDIPN(I)+1 —> IDIPN(TI) 


IDIPN(I) +> Jd 












YES XHT(I,dJ) —> XH(T) 


YHT(I,J) —> YH(T) 


ntze $0 


NO 


COMPUTE 
ERB and ERT 


TIMET+HTE(I,J) 


+TDD — TIMET 





HBRG(I,J)+ERB ~> Tl 
HTE(I,J)+ERT — 72 





XH(I)+T2ucos(Tl) — XH(I) 
YH(I)+T2sesin(Tl) —> YH(I) 
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HSM (D) 
ze @1l ¢ 


T2/HSPDT > Tl 
TIMET+TDD+T1l —> TIMET 


3 —> NREVT 
1 —» NRUNT 





CALL SNE 
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Subroutine HDS. Subroutine HDE (Detection Event) is called 
into action whenever a helicopter arrives at a new dip station. 
Using the information previously recorded by the submarine and 
helicopter maneuver events, the relative positions of helicopter 
and submarine are determined at the beginning of each sonar sweep. 
On the basis of this target range and helicopter detection range, 
a determination is made as to whether the submarine has been de- 
tected. 

Additional functions performed by Subroutine HDE are as 
follows: 


(1) Each time a helicopter dips the dip is recorded. When 
all helicopters have completed their assigned number of 
dips, or the submarine is detected, a flag is set which 
Signals the end of the play. 


(2) If the play terminates with a submarine detection, the 
elapsed time from when the helicopters arrived at datum 
until detection time, and the number of the detecting 
helicopter is recorded. This information is used by 
the main program in computing detection statistics. 


(3) A record of the time of the next intended submarine 
movement is maintained. If a submarine maneuver is 
to take place during any helicopter dip, the detection 
event is discontinued after the sonar sweep during 
Which this maneuver is to occur. After the submarine 
maneuver has been executed the helicopter dip is re- 
Sumed at the same game time at which the interruption 
occurred. 


(4) If the probability of contacting a non-submarine target 
ig greater than zero a random number is generated and 
compared with this probability. On the basis of this 
comparison and the maximum time to pursue non-sub- 
marine contacts input to the program, appropriate de- 
lays are computed and added to dip time. 


(S) On all but the last dip for each helicopter, a maneuver 
event is stored for the helicopter currently dipping. 
The time the maneuver is to take place is determined 
by adding dip time to the game time at which the dip 
commenced. Dip time is the sum of all sweep times, 
non-submarine contact delay time and the time required 
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for the helicopter to raise the transducer and make 
the transition from hovering to forward flight. 
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Subroutine HDE - Detection =vent (General) 


DETERMINE WHICH 


HELO INVOLVED 





IS Seas TeLO 


COMPLETING A av INCREMENT 

PREVIOUSLY 

INTERRUPTED DIP DIP COUNTE 
YES 


DO ONCE FOR EACH 
(aa) SONAR SWEEP 





COMPUTE TARGET 
POSITION AND 
RANGE TO TARGET 
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He (G) 











2101L-e 
ADD SWEEP TIME 
TO GAME TIME 
DID SUBMARINE 
MANEUVER DURING NO 


THLS SWEEP 






TARGET RANGE 
> MINIMUM 
DETECTION RANGE 





YES NO 











RECORD DETECTION 
AND TIME THE 
DETECTION OCCURRE 


SNE 
DETECTION EVENT 
FOR THIS HELO 
AT CURRENT 
GAME TIME 






EXIT 











SET THIS HELO'S 
FIRST DIP COUNTER 
TO NUMBER OF 
THIS DIP PLUS ONE 


LO3 





IS THIS THE NO 
LAST SWEEP 


YES 


IS THIS THE LAST \NO 
DIP FOR THIS HELO 


YES 


AVE ALL HELOS 


COMPLETED 


LAST I 


YES 


RECORD END 


OF THIS 





CALT 


PLAY 


NO 
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IS PROBaBILITY 
OF NON-SUBM ARING NO 
CONTACT > ZERO 










DID HELO GAIN 
NON-SUBM ARINE 
CONT AG 


COMPUTE ADDITIONAL 
DIP TIME REQUIRED 
FOR PROSECUTION OF 

NON-SUBMARINE CONTACT 








SNE 
MANEUVER EVENT FOR 
THIS HELO TO OCCUR 

AT GAME TIME PLUS 
TIM? TO RAISE DOME 








Subroutine HDE - Detection Event (Detailed) 


ENTER OD 





\ (XH(I)-xT)°+#(YH(1)-yr)* 
=< RTT 





NO 


RTT * DRNG(I,d) 


YES 






ree. <= Tir NO 
=~+{257 > SHORTR 
YES 


TIRST =@ Tl 
r= TIM» = Te 














1 -—>» NOVER 
NDET+1 —> NDET 
Oo + LDIP 


NSWP(I) > J TieeSWT(J)-DAT Se Tl 







XS+VXS%T2 <> XT 
YS#VYSe#Te —> YT 


= 
Nd ty 
Wt 

eee” 





HO J>yres 4 Tel =. Ff 





YES 
J 3 NRSW = DIPN(I) > nrprps)He <0, 
10 YES 






J+l —> ISHP(I) 






Tl —> TIMET 
5, => NREVI 
CALL SNE 


1 => NOVER 





RETURN 
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Hoe (D) 
4 Lr 3 







RX => POT 


FOT 4 FOPR T1+TRD —> TIMET 


YES 


(FCTs# FCO TMAX)/FOPR => FOT 2 —> NREV 
T1+TRD+FCT —> TIMET CALL SNE 
1 —> NSWP(I) 





LOT 





Subroutine PRINT. This subroutine comprises the Input-— 
Output section of the program. It is called at the beginning 
of each run and at the completion of the last run. Input data 
cards are arranged in groups headed by a card having a Hollereth 
field name punched in the first six columns. This name signals 
the information to be read from the succeeding cards. This sys- 
tem allows the user to change any number of input quantities be- 
tween runs without preparing data cards for the complete set of 
inputs. After the inputs for each run have been read into the 
machine, the input data are printed out and output data headings 
are printed. 

u ata and rules for da card pre ation. Three types 
of fields are used; bxternal Fixed Point (F), Integer (I), and 
(A). All (F) fields are F1lO.5 and when this field is specified 
data may be punched using either one of two methods: (1) with- 
out the decimal point, requiring that the last digit prior to 
the decimal point be punched in the fifth column of the field, or 
(2) with the decimal point, in which case the number may fall any 
place within the field. 

EXAMPLE: To read in the numbers 1023.65 and 20.658 using 
&a 2F10.5 field: 
(1) Without the decimal point: 
6 


Column 1 Fe a 15 li 






(2) With the decimal point: 





When an (I) field is specified the last digit of the number must 


fall in the last column of the field. The (A) field is used for 
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reading in literal characters and the letters or numbers to be 
read may fall any place within the six columns of the field. 

All arrays within a group are read in an alternating se- 
quence, i.e., If SWT (I) and PD(I) are to be read from a 
6F10.5 field, I = 1,4 the following order is used. 


Card Columns Name 
1 1-10 SWT(1) 
11-20 PD(1) 
21-30 SWT(2) 
31-40 PD(2) 
41-50 SWT(3) 
51-60 PD(3) 
e 1-10 Swr(4) 
11-20 PD(4) 


Two dimensional arrays (Array (I,J)) are listed on the card in 
order of increasing J and a new card is started for each increase 
in I. If more than one array is read in a group, the names are 
alternated as for one dimensional tables. In the following in- 
structions only the first elements are shown; succeeding elements 


Will be in the order shown. 


ata Card Forma 





GROUP I 
Card Columns Field Name Remarks 
1 1-6 A6 # DATE 
15-16 I10 NDAY Current date. 
25-26 MONTH Number of current month. 
35-36 NYEAR Last two digits of current 


year. 


This Group must be the first card and must appear only once in a 
series of runs. 


+t 
Literal characters. 
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GROUP II 








Card Columns Field Name Remarks 
1 L=-5 A6 *#NRS AM 
7-16 110 NRS AM Number of samples. 
26 KDET Always equal to one when 


included. Signals that 
detection statistics are 
to be printed at end of 





27-36 IRNG run. (See output data) 
2 1-6 12I6 NSSIZE(1) Number vs entries 
and addi-7-12 NSSIZE(2) must agree with NRSAM. 
tional ‘ ‘. 
cards if. ° 
required. . 
GROUP III 
Card Columns Field Name Remarks 
a 1-6 A6 *NHISPD Submarine speed and depth 
limits. 
7-16 110 NHISPD Knots 
17-26 LOSPD Knots 
27-36 MAXSD Feet 
37-46 MINSD Feet 


This group is normally not included if submarine maneuvers are 
predetermined. However, it is legal to include both groups III 
and V in the same set of data. The data to be used for any sub-e- 
sequent run is then determined by the value of NRANSS. 


GROUP IV 
Card Columns Field Name Remarks 
™ 1-6 A6 *NRANSS 
16 L6 NRANSS Equals one when included. 
22 NR ANC Equals zero if random number 


is to be addedto each sub- 
marine course, one otherwise. 


This group need not be included unless submarine maneuvers are 
predetermined. The one exception is if both Groups III and V 
have been read on the same run. In this case NRANSS and NRANC 
must be read in each time the submarine mode of operation is to 
be changed. <A zero in column 16 will result in random sub- 
marine maneuvers. 
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GROUP V 








Cag"a columns Field Nane 
fd L=-5 A6 #NRM AN 
7-16 110 NRM aN 
C 1-LO 4F10.5 SMTIME(1) 
11-20 DEPTH(1, 
21-30 SCUS(1) 
41-40 SSPD(1) 
5) 1-10 SMTIME(2) 
and ag . : 
required ° 
GROUP VI 
Card Columns Field Name 
ae 1-4 A6 #NRSW 
7-16 T10 NRSW 
C L-/,0 6F10.5 SWT( 1) 
11-20 PD(1) 
21-30 SWT(2) 
31-40 PD(2) 
41-50 SWT(3) 
51-60 PD(3) 
GROUP VII 
Card Columns Field Name 
pr 1-3 A6 * TDD 
2 1-10 2F10.5 TDD 
11-20 TRD 
21-30 FC PR 
31-40 FO TM AX 
GROUP VIII 
Card Columns Field Name 
1 1-4 Ab #NRHS 
7-16 110 NRHS 
17-26 NRDIPS 
27-36 HSPD 


ei 


Remarks 


Number of submarine maneu- 
vers. 

Minutes (Time the maneuver 
is to take place) SMTIME(1) 
must always be 0.0. 

Depth. (feet) 

Course. (degrees) 

Speed. (knots) 


Remarks 


Number of sonar sweeps 


Minutes to complete sweep. 
Transducer depth. (feet) 


Remarks 


Minutes to lower transducer. 
Minutes to raise transducer. 
O= MOPR 2 170 

Maximum time to classify 
non-submarine targets. 
(Minutes). 


Remarks 


Number of helicopters. 
Number of dips for each 
helicopter. 

Helicopter speed (knots). 





2 


V-10 6FLO.5 


and addi- 
tionel cards 
as needed. 


New 
Card 


Card 
1 


L1L-20 HTE(1,1) 
21-30 HBRG(1,2) 
31-40 HTE(1,2) 
41-50 HBRG(1, 3) 
1-10 6F10,5 HRRG(2,1) 
11-20 HTE(2,1) 
GROUP IX 
Columns Field Name 
1-4 AG *¥XDAT 
7-16 I10 XDAT 
17-26 YDAT 


EBRG(1,1) 


Bearing of iivrst Heli- 
copter's first dip sta- 
tion. (Degrees relative 
to 000). 

Jump distance (Yards) 
Degrees relative to 
HBRG(1,1) 


Degrees relative to 
HBRG(1,2) 


Remarks 


X and Y coordinates 
of datum. 


This group may be omitted in which case the datum coordinates 
will be (0,0). 


Card 
i. 


Card 





a 
2 


GROUP X 
Columns Field Name 
L=5 A6 * SON AR 
7-16 I10 NS 
1-10 6F10.5 FOMMU 
11<20 FOMVAR 
21-30 DL 
31-40 ‘hy 
41-50 FB 
51-60 TS 
GROUP XI 
Columns Field Name 
1-3 A6 * DAT 
1-10 F1O.5 DAT 
11-16 A6 N AME 


IL 


Remarks 


sea State 


Average sonar figure of 
merit. (Decibels) 
Variance of FOM. (db) 
Layer depth. (Feet) 
Temperature in the 
layer. (Fahrenheit) 
Sonar frequency. (kc.) 
Target strength. (db.) 


Remarks 


Time late. (minutes) 
Alpha numeric characters 
(Name of search plan) 





This group must appear in every run and must not be placed 
ahead of Groups I or II. 








GROUP xII 
Card Columns Field Name Remarks 
1 1-2 AO *GO Signals program to print 
out inputs and commence 
new run. Must appear at 
the end of each set of 
data. 
GROUP XIII 
Card Columns Field Name Remarks 
Z 1-5 A6 *#FINIS Signals program to compute 


detection statistics at the 
end of the last run and 
stop the machine. Need 
not appear if detection 
statistics are not desired. 
Output data. The output statistics reflect the writer's 
experience in ASW Helicopter squadrons. No assurance is given 
that they are suited to any other potential users requirements. 
It is relatively easy, however, to write the coding necessary 
to compute and print any desired output. 
The following output data is printed at the completion of 
every sample: 


(1) Sample size. 


(2) Nurber of detections. This is the primary measure 
of effectiveness associated with this simulation. 


(3) Percentage of plays resulting in detections. This 
is included principally as an aid in comparing 
samples of different size. 


.. ~ Number of detections . 
Percentage = Sauter ef fe 100 (6.5) 


(4) Average time to detect when a detection occurs. 
The time to complete those plays which do not end 
with a submarine detection are not included in this 


ded 3 





where: 


average. Time to detect is measured from the time 
the helicopters arrive at datum. aA problem en- 
countered in fleet ASW operations is when to aband- 
on a search once begun. Average time to detect 
should give the user some indication of the time 

at which further search with a particular search 
plan is apt to be unproductive. 


Sample variance of time to detect. 
n 
s* = (1/(n-1)) i (Ke (6.6) 
~ 421 4 . 


= number of detections 


time to detect for the {th 


detection. 


average time to detect. 


sample variance of time to detect. 
Minimum time from datum time to detection 


Maximum time from datum time to detection. Neither 
(6) nor (7) include plays which do not end in a 
submarine detection. 


Number of detections by each helicopter. In general 
considerable discretion must be used in drawing in-= 
ferences from individual unit performance. ‘Treated 
With due caution, however, some meaningful informa- 
tion may be gained from this knowledge. This parte- 
icularly applies to asymmetrical or random station 
search plans. 


Printing of the above information is predetermined by 


the computer coding and the user has no control over whether 


1 
the data is printed or not. The succeeding statistics are 


computed and written out after the last run and at the com- 


pletion of any run determined by the user. 


(1) 


Average number of detections per sample. For example: 
If this is printed for the first time at the end of 
run two, and runs one and two each consist of two 
samples. . 


D = (2) 44 D, (6.7) 
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where:D = average number of detections per sample. 


= the number of detections for the 1 sample. 


oO 
-_ 


(2) Sample variance of number of detections per sample. 
Referring to the above example: 


m 
2 =\c 
s? = (1/3)).4-, (D, - B) (6.8) 


where: S, = sample variance and D, and D are as defined above. 


(3) Maximum number of detections taken over all samples. 
(4) Minimum number of detections taken over all samples. 
The above quantities have meaning only if all samples 
are of the same size. It is therefore, anticipated that the 
user would want to compute and print this data whenever sample 
size is changed. Note that when the data is called for on any 
run, the average is taken over all runs since the last time 


these quantities were printed. 
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Subroutine PRINT General 


ENTER 
IS THIS FIRST RUN 


Q « 
YES 
READ DATE 

| 





COMPUTE 





DETECTION 
STATISTICS 
catniiiessitioni umd | TITLE 
and 
DATE 
DETEC TION 
ARE DETECTION\ YES 
es eee STATISTICS GO TOS 
TO BE COMPUTED 


NO 


. NO .( READ NUMBER 
LAST RUN COMPLETE OF SAMPLES & 
YES __ SAMPLE SIZES 

STOP 











SUBMARINE COURSE 
SPEED AND DEPTH 
HOSEN RANDOMLY 






READ NUMBER OF | ‘© 

SUBMARINE MANEUVERS, 
TIME, DEPTH, COURSE AND 
SPEED OF EACH MANUVER 








READ HIGH AND 
LOW LIMITS OF 
SUBMARINE SPEED 
AND DEPTH 













READ RANDOM 
SUBMARINE AND 
RANDOM COURSE 
INDIC ATORS 


aleb 





PRINT (G) 2 of 5 


: 


READ NUMBER OF SONAR 
SWEEPS, DEPTH AND TIME 
TO COMPLETE EACH SWEEP 


READ TIMES TO 
ESTABLISH HOVER AND 
RETRIEVE TRANSDUCER 








READ NUMBER OF 
HELOS, NUMBER OF DIPS, 
SPEED, JUMP BEARINGS 

AND JUMP DISTANCES 


a a 





CONVERT HELICOPTER | 

| JUMP DISTANCES TO 
' TIME ENROUTE AND | 
| 


| BEARINGS TO RADIANS 


om 


‘~ READ DATUM 
COORDINATES _ 


ee i ee 


READ SONAR INFORMATION 

SEA STATE, AVERAGE FOM, 
FOM VARIANCE, LAYER DEPTH, 
TEMPERATURE IN LAYER, SONAR 
FREQUENCY AND TARGET STRENGTH 







/ 
READ DATUM 
TIME AND SEARCH 
PLAN NAME 









Lil 





NEXT CARD 
FINIS 


NO 








ERROR 
PRIGT 
OUT 


NO 
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PRIff (G@) 5 of @ 


HEADING 
FOR 
INPUT 


NEXT CARD 
LABELED GO 


YES 






HEA ING 
FOR 
OU TPU 


INCREASE RUN 
NUMBER COUNTER 


RETURN 


_ Se) Se = 





Subroutine PRIN? (Detailed) 






5 
NAME,N1,N2,N3,N4,N5 K Al 
(A6,5110) 
Nl > NDAY | \;< 
N2 —> MONTH NAME = DATE 
N3 —> NYEAR 
NO 
NO 
= Pi = N rg. KDET = 0 
1 >NRU 






YES 


Nl —> NRSAM 
(21) WNC <> KDET 


NO 






3 
NRU NSSIZE(I), 
Input Data | I = 1,NRSAM 





O —> KDET 





aS 


dal 





6 

J1l,J2 

Detection Data 
Title 

a 


NRU =~ Jl 
O -> KDET 


J? — 72 
BARD/T2 —> BARD 

( (VARD-T24BARDs+#2 ) 
/(T2-1.)) —> VARD 












6 
BARD, VARD, 
MIND,MAXD 
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Jl 
Detectiun Data 
Th tle, 






—> BARD 
=p VARD 
-> MIND 
~~ MAXD 
> JT 


PRINT (D) 


of 


G 











of 9 








PRANT (D7) 
5 Of 9 


Ni -—>NHISPD 
YES (Ne ws LOSPD 
WAMéE = NHISPD N3 ->MAXSD 


N4 =<MINSD 


NO 


YES 
Nl -> NRANSS NAME = NRANSS GO 10 
NO 
YES 
<> NAME = NRM/ Nl ->NRM AN 


NO 





5 
SMTIME(I),DEPTH(I), 
SCUS(I),SSPD(I), 











NL -—> NRSW 


3 
Sut(t).PD(1), 
L = J¥_,NRSW 





el. 





PRT (Dp) 
4 of 9 


NAME={NRES a 





NO 
Nl —> XDAT YES Ty 
NO —> YDAT —_ 
NO 
YES 2 
GO TO NAME=SONAR FOMMU, FOMVAR, 
Al “= Bb, T,F-TS 
NO 6F10.5 


Nl — NCAV " NAME=NC AV 


: a> 








> 
NC AVS(I) ,NCAVD(I) 
£=1,80a¥ 


=. ae 









NCAVS(NCAV) —>» NHISPD 
NCAVS(1) —» LOSPD 

NCAVD(NCAV) —> MAXSD 
NCAVD(1) —» MINSD 


Ze 





Pee . (Dd) 
5 OF 9 


5 
HBRG(I,J),HTE(I,d) 
J = 1,NRDIPS 


(90.-HBRG(I,1))/57.29578 
— HBRG(I,1) 


HBRG(I,d-1)-(HBRG(I,d)/57.29578) 
— HBRG(I,J) 








J 2 NRDIPS > J 


YES 


~~ 1 2 NRHS 


YES 


L23 





PRINT (D) 
6 eo. © 






-—> HSPDT 


HTE(I,J)/HSPDT 
—> HTE(I,J) 
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PRIN (D) 


RP " “Ot 9 





2 
NAME=DAT YES | DAT,NAME 


O 


G0. 70 TES (y aME=GO wRU $1 «~)}SES 
f 





Ye 





PR. (D) 











S or 9 
HSPD, TDD, FOMMU, TRD, 
FOMVAR, PRFC, FOTMAX, 
NRDIPS,F,NRSW,NS,DL,T 
(I,PD(I),SWt(I),I=1,NRSW) 
NRANSS £0 Ncoav £0 
NO 
YES /wranc S$ 0 
i.) 
NO 








Title 
(SMTIME(I),SCUS(I),SSPD(I), 
DEPTH(I),I = 1,NRMAN} 
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NRU+1l —> NRU 
Jo+l —> Jo 





RETURN 


lag 


PRINT (D) 
9 of 9 





Subroutine COMPRG. This subroutine controls the compu- 
tation of sonar detection ranges for each helicopter. The 
detection ranges are computed by an iterative process using 
the applicable zone or diffraction region equations. MThese 
regions are determined by the sonar conditions input to the 
program, submarine depth, transducer depth and the helicopter 
figure of merit. 

COMPRG is a modification of a program written by 
R. L. Klinkmer of the Applied Physics Laboratory, Johns Hopkins 
University. Subroutines SCR(KP) and SETCOF together with sub- 
routine COMPRG accomplish the sonar detection range computa- 
tions. The first two of these subroutines are used as they 


Were written by Mr. Klinkner. 
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Subroutine COMPRG (General ) 


DO FOR BACH 
(a) SONAR SWEEP 


COMPUTE COEFFICIENTS 
FOR PROPAGATION 


LOSS EQUATIONS 









DO FOR BAO 
ELICOPTER 
<8, nie ANY LAYER 
YES 
COMPUTE RANGE 
YES USING ZONE 1 


AVAILABLE PROPAGATION 
LOSS FALL IN ZONE i 
Ge 


EQUATION 


NO 










ARE BOTH THE 
PROJECTOR AND 
TARGET IN THE 
LAYER 
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THE ABOVE RANGE (IF APPLIC ABLE) 


COMPRG (G) 


(2) e or 


COMPUTE RANGE 
USING ZONE II 
EQU ATION 








AVAILABLE PROPAGATION « 
LOSS FALL IN ZONE II 












COMPUTE RANGE 
USING ZONE III 
BQU ATION 







COMPUTE PROPAGATION LOSS FOR 





USING PROPER DIFFRACTION 
REGION EQUATION BASED ON 
CRITICAL PROPAGATION LOSS 

















COMPUTE RANGE USING 
AVAILABLE PROPAGATION 
LOSS AND PROPER 
DIFFRACTION REGION 
EQU ATION 







IS THIS PROPAGATIO 
LOSS LESS THAN 
THE AVAILABLE 

PROPAGATION LOSS 














DETECTION RANGE 


IS LAST RANGE 
COMPUTED ABOVE 
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Subroutine NORM. The function of Subroutine NORM is to 
compute the value of a Normal random variable, employing the 
classical central limit theorem. Given a sequence of inde- 
pendent, identically distributed random variables Xy with 
finite mean E(X) and standard deviationO(X), then the se- 


quence Y, defined by 


Y = (X, + Xn too + Xn) - nE(X) 
ee ee ie 


Va Om) 


(6.9) 


converges in distribution to a random variable which is norm- 
ally distributed with mean Zero and variance one. [27] 

In actual practice the value of a pseudo-normally dis- 
tributed random variable is computed since pseudo-random 
numbers generated by RNG are used for the sequence X,° Under 


the assumption that the xX are uniformly distributed in the 


1 
interval [o,2] then £E(X) = # and 


iL 
Xx) =- = 
wo te 


Taking n as twelve we compute 
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(x, +e o's Xy 5) - ‘| (6akO) 


Which is approximately normal with mean zero and variance one. 
To obtain the value of a normal random variable Z, with mean L{ 
and variance g * We use the transformation 
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The following calling sequence 1s used for Subroutine NORM: 
[L—-> XMU 


Oo “= $162 
CALL NORM 
The resulting normal variate will be placed into the var- 


lable YNORM by the subroutine. 
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Subroutine NORM (General ) 
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Subroutine NORM (Detailed) 


O. =—> YNORM 


0 ALL 
RAND 
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RETURN 
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